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Abstract 

We consider branes in refined topological strings. We argue that their wave-functions 
satisfy a Schrodinger equation depending on multiple times and prove this in the case 
where the topological string has a dual matrix model description. Furthermore, in 
the limit where one of the equivariant rotations approaches zero, the brane partition 
function satisfies a time-independent Schrodinger equation. We use this observation, 
as well as the back reaction of the brane on the closed string geometry, to offer an 
explanation of the connection between integrable systems and M = 2 gauge systems 
in four dimensions observed by Nekrasov and Shatashvili. 
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1 Introduction 

The study of topological strings on Calabi-Yau threefolds has led to important insights 
in various aspects of string theory and supersymmetric gauge theories. For instance, 
the properties of topological strings on compact Calabi-Yau geometries has impacted 
our understanding of BPS states of four- dimensional charged black holes. On the other 
hand, the study of topological strings in the background of non-compact 'local' Calabi- 
Yau geometries has enriched our understanding of M = 2 and M = 1 supersymmetric 
gauge theories in four dimensions. 

More recently, inspired by the work of Nekrasov on instanton partition functions 
of M = 2 gauge theories in four and five dimensions P], it became clear that there 
should exist a refinement of topological strings for the case of non-compact Calabi- 
Yau threefolds. This is because these gauge theories can also be engineered using 
non-compact, toric Calabi-Yau threefolds, for which the topological string partition 
function is a special limit of the instanton partition function, i.e., 

Ztopids) = Z Nek (ei, e 2 ) \- ei=€2= g a , 

where e a denote the two equivariant rotations of space-time (taken to be C 2 ), and g s is 
the coupling constant of topological strings. This connection implies that there should 
be a refinement of the topological string partition function, where the string coupling 
constant is split into two independent parameters. 

Despite of the above hint, a world-sheet description of such a refined topological 
string is still lacking. Nevertheless, the above observation does motivate a definition for 
the refined A-model topological string in terms of its lift to M-theory and the degen- 
eracy of BPS states in five dimensions. In particular, the standard topological string 
partition function captures the BPS degeneracy of M2-branes wrapped over 2-cycles. 
Furthermore, the topological string coupling constant is related to the chemical poten- 
tial for the SU(2)l C SU(2)l x SU(2)r = SO(4) rotation group in five dimensions. 
The refinement of the topological string captures in addition the SU(2)r rotation quan- 
tum numbers of the BPS states, and thus the full rotational quantum numbers of the 
BPS states. 

Independently, in the context of B-model topological strings that are dual to matrix 
models, an alternative definition was proposed in [2], where the refinement involves the 
so-called '/3-ensemble' of matrix models. This is just the ordinary matrix model, but 
now with the power of the Vandermonde determinant of the eigenvalues raised to 2/3. 
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The parameter (3 is related to the equivariant rotation parameters via (5 = — ei/e2 and 
equals to 1 in the case of the usual topological string. 

One natural question is how to actually compute the partition function of the 
refined topological string. There are various possibilities. For the A-model, it has been 
proposed that the refined topological vertex, which is the one-parameter refinement of 
the standard topological vertex, computes the corresponding degeneracies [3]. Another 
approach, in the context of B-model, is to use the standard matrix model techniques 
but adapted to the /3-ensemble [2] . Finally, yet another approach is to use the fact that 
the refined topological string partition function is still expected to be a wave-function, 
and thus satisfies the holomorphic anomaly equation. In this way the computation of 
the partition function reduces to fixing the holomorphic ambiguities [U [6]. Even 
though these different approaches have not been proven to be equivalent, so far they 
give results which agree with each other and with the results of Nekrasov in the cases 
that have been checked. 

On the other hand, an interesting observation concerning the instanton partition 
function was made by Nekrasov and Shatashvili (NS) [7]. They noticed that in the limit 
where one of the equivariant parameters is sent to zero while the other is kept fixed, 
for instance t\ — > with e 2 = H fixed, a connection with certain quantum integrable 
systems emerges in the following way: The integrally stationary equation of the free 
energy 

exp(0 aj W(5;fi)) = l, (1.1) 

where 

W(a; H) = lim e\JF{a\ ei, e 2 = h) , J 7 = log Z Nek 

€l->0 

and a denotes the collection of all Coulomb parameters {a/} of the theory, gives the 
Bethe ansatz equation for the corresponding integrable system. In particular the energy 
eigenvalues of the integrable system are completely characterized by the aj satisfying 
these equations. Even though this is a striking statement, and was checked in many 
cases, an explanation of it was not offered. One aim of this paper is to shed light on 
this relation. 

Our approach is to study branes in the context of the refined topological string. It 
has been long known that the branes probe geometry in a quantum mechanical way. In 
particular, in the context of the standard, unrefined topological string, it was argued 
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in [8] that for the B-model on a local Calabi-Yau given by 

uv + H(x,p) = , 

the wave-function ty(x) of a brane whose position is labeled by a point x on the Riemann 
surface H(x,p) = classically, satisfies an operator equation 



at quantum level, with H := H(x, g s d x ). In other words, one is studying a quantum 
mechanical problem given by the Hamiltonian H. Clearly, this is the right context to 
elucidate the connection to the results of NS integrability, since the B-model mirror 
to a geometrically engineered gauge theory is indeed a geometry of the above form 
and the corresponding integrable model is specified by H defined above. However, as 
pointed out in [S] , when there are cycles in the Calabi-Yau geometry corresponding to 
magnetic charges there are further g s corrections to H itself. Therefore, in general the 
above operator relation only holds up to order g s corrections. 

In this paper we consider brane partition functions in the context of refined topolog- 
ical strings. In this case there are two types of branes, corresponding in the M-theory 
language to M5 branes wrapping a Lagrangian cycle in the Calabi-Yau and a two-plane 
z\ or Z2 in IR 4 ~ C 2 associated with the e\ and e 2 action respectively. In the B-model 
context, we show that for a /3-deformed matrix model with a polynomial potential, 
the equation (11. 2p is generalized to an equation which is exact and of the form of a 
multi-time dependent version of the Schrodinger equation, i.e., 



where fi(t) are some functions of the 'time' variables tj and the momentum operator is 
given by either p = e\d x or p = e^d x , depending on the type of brane under considera- 
tion. In the refined topological B-model dual to this matrix model, the tj are identified 
with certain moduli of the corresponding Calabi-Yau geometry. This makes concrete 
and elucidates the proposal in [7] regarding the existence of the t-deformation of the 
quantum integrable system. In the NS limit we have eie 2 — > and the time dependence 
vanishes, and we simply obtain the time-independent Schrodinger equation (11. 2p for 
the e 2 -type brane. 

Note however, even though we have found solutions to the time-independent 
Schrodinger equation in the NS limit, there is no guarantee that the wave- function 



HV(x) = 0, 



(1.2) 
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has no monodromy. Some extra quantization conditions are required in order for the 
wave-functions to be well-behaved under monodromy. On the other hand, it was known 
from [8] that taking branes around the cycles of a Calabi-Yau shifts the dual periods in 
units of g s in the usual topological string. In the NS limit we have g s = (— e^) 1 ^ 2 — >■ 0, 
and the shifts become equivalent to computing derivatives. Thus we obtain the state- 
ment that taking the brane around the cycles gives the gradient of the closed string 
partition function with respect to the closed string moduli. This gives an effective way 
to compute the closed string partition function of the refined topological string theory 
in the NS limit. 

Furthermore, to have a well-defined wave-function we need the wave-function to 
be single-valued, which by the above argument requires that the exponential of the 
gradients of the partition function have to be equal to one, leading to the relation ( II. ip . 
Therefore, we understand the observation of NS regarding the relation between gauge 
theory and quantum integrable system as a consequence of the consistency of open 
refined topological string theory. Topological string theory with branes is automatically 
the relevant system to study to understand the integrable models in the NS limit, and 
the geometric engineering of gauge theory tells us that the same physical quantities can 
also be computed using target space physics and in particular the instanton calculus, 
hence closing the circle of ideas. 

Even though the derivation we found was done in the context of matrix models, we 
conjecture that it applies to all cases where refined topological strings are defined. In 
particular, we argue that one can recover the refined closed string partition functions 
in the NS limit, simply by computing the wave-function solution to the Schrodinger 
equation and seeing what monodromies it picks up as we take it around cycles. We 
verify this general conjecture in many examples. 

The organization of this paper is as follows. In section [2] we review some earlier 
results concerning topological strings and its connection with quantum integrable struc- 
tures [8], on which this work heavily relies. This is followed by recalling the various 
definitions of refined A- and B-model topological strings in section El In section H] we 
show that for B-models dual to matrix models, the brane partition functions in the 
refined topological string satisfy a time-dependent Schrodinger equation. Taking the 
NS limit of the refined topological string, the time-dependent Schrodinger equation 
reduces to the usual time-independent Schrodinger equation as we will show in section 
|5j We then use this result to explain the observation of NS about the relation between 
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gauge theory partition functions and the Bethe-ansatz for integrable systems. Sections 
[6] and [7] are devoted to explicit examples illustrating and testing our results and con- 
jectures for the NS limit in the case of genus and genus 1 surfaces. In section [8] we 
consider B-model branes in Penner type geometries. In this case our time-dependent 
Schrodinger equation is identical to the BPZ equation for the corresponding degen- 
erate fields of Liouville, and we show that our methods provide a way to derive the 
n-point function by studying the conformal blocks with degenerate insertions. Finally, 
in section [9] we present our conclusions. 

Recently there has been an increased interest in the NS limit and several papers on 
related topics have appeared during the course of this work P [101 EH EE] • 



2 Review of Quantum Geometry in Topological Strings 

In this section we review some elementary features of topological strings on local Calabi- 
Yau threefolds of relevance to this work. Since topological strings on such backgrounds 
have been extensively studied, we will be brief. 

We will be mainly interested in the B-model topological string and we focus our 
attention on local Calabi-Yau threefolds X which are given by a hypersurface of the 
form 

uv + H{x,p) = 0, (2.1) 

where u, v £ C and x,p £ C or C*. The classical (i.e., tree level) amplitude of the 
B-model is encoded via the special geometry of X in the periods of the holomorphic 
3-form 

du , 
U! = — A dp A dx . 
u 

For this class of backgrounds, the periods of u reduce to residue integrals on u = 0. 
In detail, if we view the 3-cycles as a fibration of the circle obtained by the rotation 
(u,v) — > [ue td , ve~ td ) over a disc D, where u = on the boundary of the disc, i.e., 
dD C {H(x,p) = 0}, then the period integrals reduce to 



dpdx = pdx , 

D J 3D 

where the boundary dD is a curve on H(x,p) = 0. Thus the B-model reduces to the 
study of the periods of the 1-form 

X=pdx, (2.2) 
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along the 1-cycles of the Riemann surface 



£ : H(x,p) = 0. 



(2-3) 



Hence, we end up with a local version of special geometry which involves a 1-form on 
a Riemann surface. The pair (£, A) is often referee! to as the spectral curve. 

2.1 Quantum Mechanics and Loop Corrections 

This is the story at tree level. The natural question is what the quantum corrections 
do to this picture, i.e., what is the analog of 'quantum special geometry' ? It turns 
out that an important role is played by the notion of branes. 

Consider a B-brane given by the w-plane at u — 0, concentrated on a point x of 
E, where p = p(x) is fixed in terms of x by the condition H(x,p) = 0. The classical 
partition function ty(x) of this brane, which is a function of its moduli x, was studied 
in [13] and was found to be 



where g s is the coupling constant of the topological string. This structure is very 
reminiscent of the WKB approximation to the ground state wave- function, if we identify 
H(x,p) as the Hamiltonian of the quantum system. In fact, it was argued in [8] that 
on the B-brane phase space the variables p and x do not commute and we have the 
relation 



exactly as one has in the usual set-up of quantum mechanics. Moreover this strongly 
suggests that the quantum corrections to the partition function should make *f>(x) an 
exact wave-function for the quantum Hamiltonian H, i.e., we expect a relation of the 
form 



as the operator realization of the geometric condition H(x,p) = 0. This was proposed 
in [5] and checked in various examples. It was found that for simple cases this is exactly 
right and gives the full answer for the quantum corrections to the wave-function. For 
example, in the case of the (deformed) conifold 




(2.4) 



fa, x] = g, 



(2.5) 



H(x,p) &(x) = 



(2.6) 



H(x,p) = —p 2 + x 2 — ji , 
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the full partition function of the brane is indeed given by the corresponding energy 
eigenstate of the harmonic oscillator. 

A motivation for the above commutation relation is the following. The brane wrap- 
ping the w-plane has fields x, p living on it that capture the normal deformations of the 
brane as a function of v. The kinetic term for this brane is given by 



which leads to the fact that p and x are conjugate variables. 

A perhaps more familiar setup to motivate (I2.5P is via mirror symmetry. When X 
is a mirror to a toric Calabi-Yau manifold X*, the B-brane on the v -plane gets mapped 
to an A-brane on X* wrapping a non-compact Lagrangian three manifold, L = M 2 x S 1 . 
The theory on a Lagrangian B-brane is a U(l) Chern-Simons theory, with the classical 
action 



In the absence of world-sheet instantons, the action (I2.7P is exact. The Chern-Simons 
path integral on L computes the A-model partition function with the brane inserted 
[H]. Since L is a manifold with a T 2 boundary, the Chern-Simons path integral on L 
determines a state in the Hilbert space of the theory on T 2 x l t in a familiar way. As 
is evident from ( 12. 7p . when quantizing Chern-Simons theory on T 2 x R with M. viewed 
as time, the holonomies of the gauge field around the (1,0) and (0,1) cycles of T 2 
become canonically conjugate to each other. In fact, mirror symmetry directly relates 
the B-model variables x and p to the complexified holonomies 



In general, the world-sheet instantons will correct the classical action, but not in a way 
that affects this observation. (Namely, they shift Sqs to Sqs + S'inst, but the instanton 
generated terms are independent of time.) 

An example of this is given by the mirror of the A-model on the (trivial) background 
X* = C 2 x C*, with the S 1 in L mapped to the circle in C*. In this case, the mirror is 



where the holonomy around the S 1 maps to x by mirror symmetry. In this case, there 
are no world-sheet instantons, and (I2.6P leads to the trivial result 





(2.7) 




H(x,p) = e p -l, 



$>{x) = 1, 
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which is the correct answer for the Chern-Simons partition function on L. 

A less trivial example corresponds to the A- model on X* = C 3 . This maps to the 
local threefold geometry where the spectral curve is given by [T5J 



H(x,p) = e p + e x -1 = 0. 
The equation ( 12.61) now leads to the wave-function given by the quantum dilogarithm 

oo 

= J] (1 - eV) , 

71=0 

where we defined q := e 9s , as the solution capturing higher loop corrections for the 
open topological string. This is mirror to the A-model partition function, when one 
includes instanton corrections along with the Chern-Simons functional [TB| [T5]. 

It was found in [8] that this simple picture is also true for all the B-model geome- 
tries where E is a genus zero Riemann surface. In terms of the mirror A-model, this 
corresponds to geometries with no closed 4-cycles. However, it was also found in [S] 
that the picture is more complicated if £ is a higher genus Riemann surface. Namely, in 
general, ty(x) is a zero eigenstate of an Hamiltonian H'(x,p) which is not the classical 
Hamiltonian H(x,p) but receives g s corrections. In other words we have a quantum 
Hamiltonian ground state 

H'(x,p)^(x) = 0, where H'(x,p) = H(x,p) + 0(g s ) . 

We do expect such quantum corrections, since in general there are already normal or- 
dering ambiguities in replacing the classical variables p, x with the quantum operators. 
This shifts the question of the utility of this approach to a better understanding of 
how to compute corrections to H, which was not addressed in [8], but will be answered 
later in this paper. 

2.2 Branes and Period Shifts 

Besides the calculation of quantum corrections, there is a further issue to address. 
So far we have only discussed the partition function of the topological string in the 

1 Note that the fact that only terms ~ e nx+mp w j^ n n ^ m being integers appear in the equations for 
B-modcl mirrors is related to the fact that x and p are variables mirror to holonomies on T 2 , which 
are clearly periodic. 
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presence of a brane. How can we obtain the partition function of the closed topological 
string in the absence of the brane? 

These two issues are in fact closely related. It turns out that the insertion of branes 
affects the closed string moduli [Hj . The meromorphic 1-form ( 12.2ft acquires a first order 
pole with residue g s at the position xq where a brane was inserted, i.e., 

A ~ — — — dx + . . . . 

X — Xq 

This implies that, if we measure the change in the period of A around the position of 
the brane, we find that we pick up a period 

<f\ = g a . (2.8) 

Jx 

The content of this statement is that topological branes source the topological gravity 
fields, i.e., they change the background geometry. This interpretation is in particular 
supported, in the mirror setup, by the large N duality of the topological A-model on 
the (singular) conifold with the resolved conifold T*S 3 [T7j. Namely, putting N branes 
on S 3 changes the Kahler class surrounding the S 3 by Ng s . 

In [8] this statement was used to provide a quantum definition of D-branes in the 
closed string field theory of the B-model. Namely, fluctuations of the one-form A 
corresponding to quantum fluctuations of complex structures on X are captured by 
the so-called Kodaira-Spencer scalar field on £ as 

5X = d<f) . 

Since is a closed string field, its two point functions are proportional to gl, i.e., 

4>(x)(j)(x') ~ g 2 s ln(x — x') + . . . . 

Moreover, deformation of the geometry in (12.81) is consistent with the brane insertion 
operator ip(x) being a fermion 

related to by bosonization. This is also consistent with the period 

(b d<p i/j(x ) = g s ip(xo) ■ 

J X() 

Furthermore, if we define 4>(x) so that it includes the classical piece corresponding 
to f x A = f x p dx, then naturally (12 .4p is the classical piece of the insertion of the 
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operator ijj(x) at the point x on the Riemann surface. The fact that the B-brane is a 
fermion of the Kodaira- Spencer scalar field <j)(x) has been used in [8] to compute the 
exact B-model amplitudes on some local Calabi-Yau geometries. 

The change of moduli due to the insertion of branes means that if we consider 
a brane/anti-brane pair, and take one of the branes around a 1-cycle 7 on E before 
bringing it back to annihilate with the other brane, we have changed in the process the 
periods of A along any other 1-cycle a by the amount 



where (0,7) denotes the intersection product of the two 1-cycles. 

In particular, if we denote by a/s the periods of the A-cycles in the Riemann surface 
E and decompose the cycle 7 into a combination of symplectically paired B-cycles, i.e., 



where a denotes the collection of the Kahler moduli. This implies that once we know the 
partition function of topological strings in the presence of branes, we can in principle 
find the closed string partition function by its variation with respect to the closed string 
moduli, at least for shifts of a/ by integral multiples of g s . 

Interesting as these statements are, they leave a gap in their applicability. This 
is because for general cases, we cannot compute the brane wave-function since 
we do not know how to compute the quantum corrections to the quantum Hamil- 
tonian H(x,p). It turns out that these problems get remedied when we consider a 
one-parameter deformation of topological strings, inspired by the related deformation 
of gauge theory partition functions [I]. 

3 Refined Topological Strings 

There is an interesting one-parameter extension of topological strings in the case of 
non-compact Calabi-Yau backgrounds. Denoting this parameter by /3, the extension 
can be combined with the usual string coupling constant g s into two parameters e±, e<i 
defined by 




7 = (A I ,B J ) = 6 J I 



we find that the closed string partition function is shifted by 
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The standard topological string corresponds to the special case (3 = 1, where we have 

- ei = e 2 = g s . (3.2) 

However, unlike the usual topological string theory, we do not yet have a world-sheet 
description of the refined topological string with generic extra parameter /3. Neverthe- 
less, in various cases we do have a definition in terms of space-time physics. Let us 
now review the different definitions. 

In the topological A-model we can give a general definition for this one-parameter 
deformation in terms of the partition function of target space physics. Assuming mirror 
symmetry is general, this would be a working definition for the topological B-model 
as well. The extension in the A-model case is accomplished by lifting to M-theory. 
The addition of an extra dimension closely parallels the construction of Nekrasov of 
^-backgrounds in M = 2 gauge theories. In particular, when the Calabi-Yau geometry 
engineers a four- dimensional gauge theory, the computation of the target space par- 
tition function can be captured by the instanton calculus of [T]. More generally, the 
one-parameter extension of the topological string has a space-time analogue which is 
a global M = 2 supersymmetric system in four dimensions. 

In the next subsection we will first review the definition of A-model topological 
strings in terms of M-theory and show how this definition can be extended to the cases 
where the Calabi-Yau is non-compact. Then we will relate this to the spectrum of BPS 
states in five dimensions. We then connect this picture to Nekrasov's definition in the 
special case where the Calabi-Yau engineers a five- dimensional gauge theory and briefly 
comment on this extension for a general global M = 2 system. Finally, we discuss a 
direct B-model definition of this extension when the topological string can be realized 
in terms of matrix models. 

3.1 The A-model, its M-Theory Lift and Deformation 

The partition function of the topological A-model on a Calabi-Yau X is equal to the 
partition function of M-theory on the space [IS] 

IxS'x TN , 

where TN is the Taub-NUT space, viewed in complex coordinates Z\,Z2 ■ Furthermore, 
the Taub-NUT space is twisted along the S 1 , in the sense that going around the circle 
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rotates the coordinates z±,z 2 by 



zi -> e l6l Zx, z 2 ->■ e l62 z 2 . 

In order to preserve supersymmetry we need ei + e 2 = , and we find (13. 2p . The refined 
topological string is obtained by relaxing this constraint on the e a . However, to do so 
while preserving supersymmetry, we need an extra U(1)r symmetry acting on X. This 
can only be accomplished with X being non-compact, a condition we will assume from 
now on. 

This partition function is not easy to compute, but can be related to the spectrum 
of BPS particles in five dimensions [191 120] • Namely, the corresponding M-theory 
partition function simply computes the BPS indices 

N? ■ 

3L,3R 

of BPS states of M2 branes wrapping 2-cycles of X of class d £ H 2 (X,Z*), and with 
5*0(4) = SU(2)l x SU(2)r spin quantum numbers jL,jn- Decomposing the corre- 
sponding field into modes 



^(zi,Z 2 ) = ^ a ni,n 2 zT 



~n 2 
^2 ' 



n\,n 2 



the U{1) L x U(1) R C SU(2) l x SU(2) r acts on this by 

( Zl ,z 2 ) -> (e ti9L+eR) z ll e l{ " eL+dR) z 2 ). 

The BPS partition function is defined as a trace in the Hilbert space Hbps of a gas of 
spinning M2 branes, weighted by their total spin and charges: 

A BPS - ^Hbps\- 1 ) 1r e 

where we have used (m/,,m^) to denote the spin content {j\,j R ) of the highest spin 
state in a given BPS multiplet. In this case, the Hilbert space is just the Fock space 
of a single M2 brane and hence the partition function takes the form0 

3l 3r N d 

ZbPS = J] _ e -^l{rn L +m R +\+m) e e 2 (rn L -m R +\+n 2 ) e -a-d\ 

d,3L,3R m L=-3L m R =-j R n 1> n 2 >Q 
2 Strictly speaking, this statement is true only with large three-form fluxes turned on |21j . See also 
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where q^ = e~ 4 2 , = e L 2 -2 " and a = (a l5 02, ... ) the set of Kahler parameters. 
Note that in addition to the intrinsic angular momenta of the particle in five dimensions, 
the orbital angular momenta contribute to 2l,r as well. 

In the special case of q R = 1 we effectively trace over the SU (2) R quantum numbers. 
In this limit, the BPS partition function equals the topological string partition function, 
and we have the usual relation between the topological string and M-theory amplitude. 
For the refined version of the topological string we define 

Ztop{a>\ e ij £2) = Z BP s(a; e±, e 2 ) . 

Equivalently, we could have obtained Tbps — ^ogZ BP s from a one- loop, Schwinger 
type computation which involves integrating out M2 brane particles in a gravi-photon 
background parameterized by e\ and e 2 that is no longer self-dual [20]. This gives J-'bps 
as 

00 00 —ka-d 
d jL,3R=0 k=l 

^ e -(ei+e2)fcjfl _|_____|_ e (ei+e2)kj R \ ^g-(-ei+e 2 )fcjL _|_ _ _ _ _|_ 

X 7TS ZTHT7 tH -a- 12. n • 

(e 2 — e 2 j(e 2 — e 2 J 

(3.3) 

For some, but not all choices of X, this will also have a gauge theory interpretation in 
terms of a five-dimensional gauge theory on a circle. 

Branes in the Refined A-model Topological String 

One can extend the above formulation to include A-branes. The usual A-type brane 
lifts in M-theory to a M5 brane wrapping the mirror Lagrangian 3-cycle L in X, and 
the C x S 1 subspace of TN x S 1 [16]. A similar consideration also applies to the refined 
version of topological strings. In this case we have two inequivalent types of branes, 
depending on which of the cigar subspaces of TN, i.e., the Z\- or the z 2 -plane, the 
M5 brane wraps. For general e Q 's, the symmetry between z\ and Z2 is broken and 
the two types of branes are no longer on equal footing. At low energies, the theory 
on the brane has M = 2 supersymmetry in three dimensions. After introducing A- 
branes, in addition to M2 brane particles wrapping closed holomorphic 2-cycles in X, 
there are also M2 branes wrapping holomorphic disks in X and ending on L. The 
branes break the local 50(4) rotation symmetry to SO(2)i x 5*0(2)2. Combining 
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the S0{2) L subgroup of it together with the 5*0(2) R-symmetry, we still have an 
SO(2) L x S0(2) R = U{1)l x U{1)r symmetry available [23]. In the presence of the 
branes we consider the BPS partition function of open and closed M2 branes, keeping 
track of their spin and the relative homology class in H2(X, L). The latter corresponds 
to the fact that the open M2 branes are charged under the world-volume (magnetic) 
gauge field on the M5 brane (see [16] for more details). In this way we can define the 
open BPS partition function, analogously to the closed one. We refer to the partition 
function in the presence of the two different types of branes, corresponding to the M5 
brane wrapping the £0=1,2 plane, as \l/ a . 

3.2 fi-Deformation of J\f = 2, d = 4 Theories 

The so-called f2-background for four- dimensional M = 2 supersymmetric theory, where 
the parameters e 1)2 correspond to the equivariant U(l) x U(l) action on the Euclidean 
space-time C 2 = M 4 , has been considered in [H [211 I2H], based on earlier work [261 
[27]. More precisely, one considers a deformation of the theory where one replaces the 
adjoint- valued scalars $ by 



corresponding to adding generators of infinitesimal rotation along the two complex 
planes of C 2 , with coordinates 2^2, accompanied by an R-symmetry twist. Since $ 
gives masses to charged fields in the theory, this deformation effectively adds a mass 
to the modes of the fields in the theory according to their transformation properties 
under the U(l) x U(l) action. The partition function of the theory in this background 
defines a refined partition function Z/v= 2 (ei, e 2 ). 

When the M = 2 theory in question is a gauge theory, the fi-background allows one 
to compute the partition function explicitly by performing integrals over the instanton 
moduli spaces. Moreover, if we consider a five-dimensional gauge theory with an ad- 
ditional S , the partition function can also be computed using the instanton calculus. 
Furthermore, this maps to our M-theory construction above in the cases where the 
local Calabi-Yau X engineers the corresponding 5d gauge system, and thereby makes 
contact with the refined A-model topological string. In fact it was explicitly verified in 
[251122] that in the special case f)3.2p . the partition function Zj^/ =2 (ei, e 2 )|- (E1 = (E2 = 9a agrees 
with the topological string partition function Z t op{g s ) on the corresponding Calabi-Yau 
manifold. Furthermore it was checked that a refined version of the topological vertex 




(3.4) 
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not only reproduces the SU{2) L x SU(2) R quantum numbers of BPS states, but it also 
agrees with the full Nekrasov partition function of the M = 2 theory [3] . 

This viewpoint on deformation also extends to the open topological string. The A- 
branes can be given a purely gauge-theoretic formulation as surface operators, studied 
recently in [301 ED E2l [33] , for example. For general deformation e a 's, the symmetry 
between Z\ and Z2 is broken and we obtain two types of branes. The breakdown of 
symmetry can already be seen at the classical level. 

From the string theory perspective, the surface operators are described by the 
branes of the theory. The world- volume theory on the brane wrapping the ^i^-plane 
starts out as an M = (2, 2) supersymmetric theory in d = 2 with a superpotential 
W(x), where x is a chiral superfield. Since the equivariant action is in essence a kind 
of Kaluza-Klein reduction that effectively gives the z a plane a volume proportional to 
l/e Q , the classical partition function of the brane wrapping the z a plane becomes 



In the case of the B- model topological string, as discussed in section the correspond- 
ing superpotential is identified with 



for the two types of branes wrapping the z^-plane. 

3.3 The /3-Ensemble as a Refinement of the Topological B-model 

The above discussion was geared towards A-model topological strings. By employing 
mirror symmetry it also gives us an answer for the B-model in principle. However, 
given the simplicity of the B-model, it is convenient to give a direct definition of the 
refined B-model. Indeed a proposal for such a definition has been put forth in [2] in 
terms of a certain deformation of matrix models. 

Large N matrix models provide an alternative description of topological B-models 
on a class of geometry and an excellent testing ground for the ideas we just reviewed. 
In view of its importance in the rest of the discussion, we will now review the relation 




a = 1,2 . 





leading to a WKB-type wave-function 




(3.5) 
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between random matrices and topological strings. We will start by first reviewing the 
usual B-model and focus on the refined case afterwards. 
For the type of geometries described by the local curve 

H(x,p) = -p 2 + W'(x) 2 (3.6) 

with W(x) a degree g + 2 polynomial, there is a conifold singularity near each critical 
point of the potential W(x) which we can blow up into a P 1 . The B-model topological 
string in this background has a matrix model description (Mj. Its partition function is 
given by 

Z = JdQe-h™*™ , (3.7) 

expanded around a given distribution Ni of eigenvalues of the matrix $ among the 
g + 1 critical points of the potential W(x). 

This matrix model can be derived from the B-model [3^| |8] as follows. If we cut 
the geometry into two halves, corresponding to writing (j3.6jl as 

-H(x,p)= (p + W'(x))(p-W'(x)) = H L (x } p)H R (x,p), 

the branes on the P 1 are obtained by gluing non-compact branes on Hlh(x,p) = 
across their boundary. From Hl, with branes inserted at x — for i — 1, . . . N we get 

^( Zl )^(z 2 )..^(z N )) L = e-^ w ^ Y[{zj-z k ). (3.8) 

j<k 

The interaction term comes from the two-point functions of ip{zi) with each other, while 
the potential term is the wave- function of a single D-brane as we explained before. Since 
this is a genus zero Riemann surface, the result is exact as we will argue later. Similarly, 
from Hr(x,p) we end up getting another copy of this, as the orientation of the branes 
is naturally opposite. Setting the values of equal on both sides and integrating them, 
we indeed obtain the partition function (13. TJ) , which reads in the eigenvalue basis 

Z= f d N zH(z i -z j ) 2 e-^ w ^. 
J i<j 

It was conjectured in [35| 131] that the geometry with the branes go through a 
transition to a smooth geometry where the P^s are replaced by S 3, s with size 

Hi = g s Ni . 
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More precisely, we now have the geometry given by 

H(x,p) = -p 2 + W'(x) 2 + f(x) , (3.9) 

where the correction (deformation) f(x) term is a degree g polynomial and determined 
by the distribution 

N = Ni + . . . + N g+1 , 

of the N branes over the g + 1 P^s given by critical points of W(x). The open-closed 
duality states that the matrix model, in the large N 't Hooft limit 

N — > oo, g s — >■ , Ng s fixed , 

gives a description of the B-model on the geometry (13. 9p after the transition. The 
curve 

£ : H(x,p) = -p 2 + W\xf + f(x) = (3.10) 

is now the spectral curve of the matrix model, from which the higher g s correction of 
the matrix model partition function can be computed [36] . 

From the same argument, we can also deduce that the insertion of a non-compact 
brane corresponds to the insertion of the characteristic determinant on the matrix 
model side and leads to the brane partition function 

Z bm ne(x) = (4>{x)) = e~Ts W{x) J d$e-Ts™^&e%{$-x). (3.11) 

The generalization of the above picture to the refined background is proposed in [2] . 
For (3 7^ 1, or e\ + e 2 7^ 0, the matrix model obtains a different measure while retaining 
the same potential. The generalization simple changes the power of the Vandermonde 
determinant to 2/3, together with a rescaling of the coupling constant in front of the 
potential W, i.e., 

Z = I d N zH( Zl - Zj re- 2 ^^ w ^. 
J %<j 

In terms of the variables e\, e 2 this translates into 

Z= j d N z H( Zl - Zj )~ 2 ^ e~^ W{zi) . (3.12) 
J i<j 

In particular, for logarithmic potentials relevant for the discussion on relations to 
J\f — 2, d — 4 gauge theories, the /3-ensemble matrix model integral takes the form of 
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a Coulomb gas representation of the conformal blocks of Liouville theory with back- 
ground charge given by 0. 

This deformation can also be understood from our space-time perspective. Namely, 
if we consider branes wrapping the z a plane, in the space time C 2 , rotated by e a , 
the classical partition function becomes (I3.5p . instead of (12. 4p . So, in deriving the 
matrix model, as we reviewed in section |2j we expect to simply replace the fermion 
ip(x) = exp(4>(x)/g s ), by the operators 

i/) a (x) = exp((j)(x)/e a ), a = 1,2. (3.13) 

On the other hand, the scalar field 4>(x) still has the same correlation functions as 
before (this is the essence of the Coulomb gas formalism), which is consistent with the 
fact that the quantum mechanics of the closed string is unchanged, as we will later 
discuss in the context of holomorphic anomaly. This means that the correlator (I3.14p 
computing the partition function of branes on the halved geometry changes to 

{M*i)M**) ■ --M*»))l = 11^ - Zj )^e-^ W{Zi) , (3.14) 

i<j 

which leads to the partition function (I3.12p upon gluing, using the fact that g 2 s = —e^. 
From the saddle point equation 

W'{ Zi ) = -e x £ -±- , 
we see that the relevant 't Hooft coupling is given by 

tM=€!N. 

We see that the amount by which the branes change the geometry around them now 
depends on the type of the brane. More generally, from (13.131) we expect that inserting 
an 61(62) brane at a point Xq on the Riemann surface deforms the geometry by 62(61). 

<f 90^(x o ) = -^(x ). (3.15) 

We can also easily describe non-compact D-branes for this matrix model. The two 
different kinds of branes 

ip* a (x) = exp(-0(x)/e a ), 
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correspond in the matrix model language to the operators 

e^ (l) • det(x - , a = 1,2, 

respectively. Such an explicit realization of branes in the matrix model allows us to 
study them very directly. We will come back to this theory in the next section, where 
we will derive the exact, time-dependent Schrodinger equation that they satisfy. 
Before we go on, note that we can rewrite the partition function (13.121) as 

Z= jd N z H(z z - 2j )2 e ^(^^(,) + 2( £l+Q )E 1<3 io g (,-, 3 )) (3 16) 

and view the change of the measure from the usual Vandermonde squared as adding a 
non-local operator with a coupling constant 

h := ei + e 2 . 

Writing the /3-ensemble partition function as in (13. 161) . makes manifest that apart from 
the usual genus expansion in g s , the free energy J 7 = logZ has another expansion in 
terms of the parameter h defined above. We end up with the double expansion |f| 

F(S;e 1 ,e 2 )= ■ (3-17) 

g>0,£>0 

The expansion (13. 1 71) and its recursion relation has been studied in [371 EHl 139] . 

3.4 Refinement, Topological Strings and Quantum Mechanics 

One striking aspect of topological string theory, is that quantum mechanics makes 
two independent, though related, appearances in the theory. On the one hand, as 
we reviewed above, the open topological string partition function is a wave-function 
with the Riemann surface as the level set of the Hamiltonian. As we will see later, 
this continues to be true after we turn on the two independent parameters ei^. On 

3 Note that the refined matrix model free energy (|3.17|1 generally possesses an expansion into even 
and odd powers of h. From the M-theory perspective reviewed in section 2, we expect that the refined 
topological string partition function actually possesses an expansion into even powers of h only. This 
simply comes from the fact that the BPS states fit into complete spin multiplets (|3.3|) . Apparently, 
the /3-ensemble breaks this symmetry. However, as anticipated in [5], the symmetry can be restored 
via appropriate redefinition of the (deformation) parameters as we will see later both in the general 
discussion and in the explicit examples we will consider. 
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the other hand, the holomorphic anomaly of [40] implies that the closed topological 
string partition function is also a wave-function. In other words, it is a state in the 
Hilbert space obtained by quantizing H 3 (X, C) [H]. This remains true for the refined 
topological string [4]. As we will explain, this fact can be understood from the way 
branes deform the geometry and the monodromy transformation of the brane partition 
function. Moreover, in section [5] we will see how the monodromy properties of the open 
topological string partition function can be used to derive differential (or difference) 
equations that the closed string partition satisfies. 

Given a symplectic basis of H 3 (X, Z) with Aj n B J = 5j , the periods of the holo- 
morphic three form u parameterize the phase space H 3 (X, C). In our case, this is 
equivalent to periods of the one-form A on the Riemann surface S 



Classically, a T D and aj are not independent, but satisfy the special geometry relation 



where we used hats to distinguish the operators from their expectation values. In 
particular, (13 . 18[) is a semiclassical approximation to this. 

The above quantum equation as well as its generalization to the refined topological 
strings can be derived by considering the way the branes deform the geometry. As 
we have seen in section 13.31 an e a -brane deforms the period of d<f> around any cycle 
surrounding it by f)3.15p . 

Consider creating an e Q -brane/anti-brane pair at a point on E, and taking one of 
the branes around a cycle 7 before annihilating with each other. Taking the branes 
around ja — Yli ^ ^1 does not change the expectation value 




A. 
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(Note that § dip is the quantum generalization of § pdx.) However, such a monodromy 
around '-/a changes the phase of the partition function by 



M 1A : Z top (a) -> exp f ^XX 0/ 1 z to V {a) , 

since this change of phase is induced by the change of phase of 

(...exp(-0(x)/e Q )) , 

as we transport the brane around. This is to be contrasted with the situation if we 
take the branes around a 7# = mjB 1 cycle 

M 1B : Z top (a) -> Z i0J) (a + — m) = exp — Z tojJ (a) , 

e a V e a ^y' da i J 

as a consequence of (I3.15p . generalizing the shift reviewed in section I2~2l Now, consider 
writing Z in the terms of the dual variables associated with the B-cycles instead. 
In this dual basis the monodromy around the cycle 7^ acts on Z as a multiplication 
operator by exp(^ 7 j-mia^), and that around the cycle becomes a shift operator. 
The consistency of the two dual pictures requires a T D to be realized as 

1 2 d 
aD = 9 ^ 

acting on the partition function in the a/-basis. We conclude that the commuta- 
tion relation ( I3.19P between the operators a/, a J D also holds for the refined topological 
string, and therefore the closed string partition function Z is indeed a wave-function 
on H 3 (X, C) for arbitrary 61,2- The only effect of the (3 deformation is to change the 
unit of the shift and hence the form of the wave-function. This is consistent with the 
observation made in [U [5] that in known cases the refined topological string partition 
function still satisfies the holomorphic anomaly equation of |40j . 

4 Matrix Models and Schrodinger Equations 

In this section, we will show that we can derive from matrix models a multi-time 
dependent Schrodinger equation for arbitrary (3, satisfied by a brane probing a Riemann 
surface (I3.10p . with 

9+2 

W(x) = J2tnX n (4.1) 

n=0 
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and f(x) a polynomial of degree g. We will sometimes take g to infinity, so that W(x) 
becomes a formal sum. The relation between quantum geometry and the /3-ensemble 
matrix model has been discussed in [121 H3j HI] in terms of the resolvent of the matrix 
model. Here, for us we the natural object to study is the brane partition function for 
which the quantum geometry of the matrix model becomes manifest. 

We will also see that the time- dependent Schrodinger equation satisfied by the 
brane wave-function can be rephrased in the form of the BPZ equation [45] satisfied by 
a correlation function in two-dimensional CFT with a degenerate operator insertion, 
similar as for surface operators in the Liouville context [9]. 

4.1 Time-Dependent Schrodinger Equation 



From section |3~3| inserting an e a -brane, the refined topological string partition function 
on this geometry becomes 

Za (x) = e^ x) j d N zH( Zi - Zj y 2 ^ l[(x - ztf^ e -£^W . (4.2) 

i<j i 

To derive the equation satisfied by (14. 2p . consider differentiating it with respect to x. 
For simplicity of the derivation, we will remove the classical piece e 7 ^ w<yX ^ > for the time 
being and restore it later. Moreover, we will denote by h = e±/e a the power of the 
determinant det(x — $) = Yli( x ~ z i)- Different values of h correspond to different 
types of branes. 

It is easy to see that differentiating is the same as inserting the following function 
inside the intergral 

d 2 v ^ h 2 x ^ h 



dx 2 (x — Zi)(x — Zj) [x — z^ 2 



^ 2h 2 h 2 -h 



. (X - Zi)(Zi - Zj) ^ (X-Zi) 



2 ' 



At the same time, from the loop equation 

= £ [d N z^- ( -^—0(z)e-^ wM Jl(z j - z k )~ 2 ^A (4.3) 
i=i J ozi^x Zi , <k J 

with the choice of operator O = Yli( x ~ z i) h i we obtain the identity 
1-h 2 f^W'izi) ^ 1 



E 




(x-Zi) 2 e 2 \' L ^x-z i ^ (x - Zi)(x - Zj) 
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under the integral sign. As we will mention in more details later, the above loop 
equation has the meaning of the Ward identity of a quantum symmetry of the matrix 
theory. In order to have a simple operator equation, we would like to cancel the terms 

7 — ^~T2 and Y]^a ? 77 7 and the above equations show that this is possible 

exactly when h takes the two values h = ei/e a corresponding to our two types of 
branes! 

Reincorporating the classical piece e 7 ^ w ^ and putting everything together, we find 
that Z a satisfies an operator equation (under the integral) of the form 

- e 2 a d 2 x + W'{xf + e a W"(x) + f(x) = , (4.4) 

with 

w ^ Zt-X 

i=l 

The effect of inserting f(x) is the same as acting by a linear differential operator 

f(x)=g 2 s J2x n d {n)l 

n=0 

with 

9+2 d 

d {n ) = Yl H *fr • 

k=n+2 Wk-n-2 

Here we differentiate the matrix model partition function with respect to the coefficients 
tk of the potential. We have also set = — 

In summary, we find that the brane partition functions Z a (x) satisfy a linear dif- 
ferential equation 

{-e 2 a d 2 x + W'(x) 2 + e a W"(x) + />)) Z a (x;t) = 0, a = 1,2. (4.5) 

We emphasize that we have not taken any limits here - the equation is exact. 

This equation is in fact a multi-time dependent Schrodinger equation, with the 
Riemann surface playing the role of a time- dependent Hamiltonian. In order to see 
this, we can proceed as follows. To begin with, Z a (x) contains both open and closed 
string contributions, since it corresponds to an unnormalized expectation value. Corre- 
spondingly, the time-dependent Schrodinger equation we got has no reference to filling 
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fractions, or the choice of the background, as they do not enter the matrix integral ex- 
plicitly. For many purposes, the normalized expectation value ^f a (x) is a more natural 
quantity 

t / \ Z n (x) 



Z 

Here Z is the matrix model partition function without the brane, in the sector corre- 
sponding to specific filling fractions 

N^(N u ...,N g+1 ), 

around which we expand in the 't Hooft expansion, Z = Z{p,j\ e±, e 2 ) where \ii = Njei 
are the 't Hooft couplings. The normalization induces explicit dependence on the 
background. It follows that ^ a {x) is the purely open string partition function and 
satisfies 

(-4^ + W\xf + f(x) + g\ £ x n d {n) ) * a (x) = , (4.6) 

n=0 

where 

9 

f(x)=J2x n b n (4.7) 

n=0 

is now a polynomial with some coefficients b n that parameterize the complex structure 
moduli, and thus implicitly contain the choice of the background. More precisely, 
writing logZ = Tdosed/g^, we have that 

b n = <9( n) 7 ' dosed + e a (n + l)(n + 2)t n . (4.8) 

Notice that we have incorporated the term proportional to W"(x) into our definition 
of the polynomial f{x). This shift has origin in the ordering ambiguity in quantization 
and corresponds to the \h shift in the ground state energy of the harmonic oscillator 
when we consider the Gaussian potential. 

To understand the meaning of this consider the planar limit 

ei — )• , Nj — > 00 with /i/ = iV/ei fixed . 

We immediately see from (14. 6 p that a particularly interesting limit is to take this planar 
limit while keeping 

e 2 = h fixed . 



26 



In this limit the time dependence of the differential equation (14. 6 j) drops out and we 
are left with an interesting time-independent Schrodinger equation for the e 2 -brane. 
We will devote section [5] to the discussion of this limit. 

On the other hand, the theory becomes classical if we in addition take e 2 to zero. 
In this limit, the Schrodinger equation becomes the Riemann surface equation (13. 9p . 



Here the coefficients b n are fixed in terms of the filling fractions, either by requiring 



around the cuts that open up from the critical points of the potential W(x), or equiva- 



and moreover, the times begin to flow, as the Schrodinger equation is time dependent. 
The b' n s started out as complicated, e a dependent function of the 't Hooft couplings, 
\ii determined from (I4.8p . However, since they parameterize the closed string moduli 
space equally well as /i/s do (since there are as many of them as there are cuts and 
the B-periods), we can use them to parameterize our ignorance of the (in general) 
complicated closed string amplitude that underlies them. Later on, we will see that 
from the solutions to the Schrodinger equation \l/(x; {b n }) we can in fact determine 
^closed, at least in the NS limit. 

4.2 Virasoro Constraints and Hidden Conformal Symmetry 

The matrix model has a conformal symmetry, which is well known in the unrefined 
case, and which survives refinement. It leads directly to the Schrodinger equation. 

Just as in the usual matrix model [IB] , in the /3-ensemble matrix model with poly- 
nomial potential W(x) we can identify the scalar field as 



p 2 + W'{xf + f(x) = 0. 




lently from the classical pre-potential J-'^(fx), by b n = d( n )J-( >. For general e 1 ,e 2 , the 
Riemann surface becomes quantum, as p becomes an operator 







P = e a —, 
ox 



<f>(x) = g^Wix) 



Tr log (a; — $) . 



Using the operator equation 



Tr$ 



e 2 d 



2 dt 
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we get the standard mode expansion of a free chiral boson 

OO OO rj 

m*) = 97 1 E kt ^ k ~ l + % E x ~w ■ 



2 ^ dt k 

k=l k=0 



The symmetry of the matrix model corresponding to re-parametrizing the eigen- 
values hides a conformal symmetry. Concretely, the Ward identity 



N 

° = E/ rf 



dzi V x - Zi 



ensures the invariance of physical quantities under reparametrizing the eigenvalues via 
arbitrary polynomial functions. 

Let T(x) be the energy-momentum tensor 

T(x) = =: d(j){x)d(t){x) : «9 2 0(x) , 

Qs 

k 38 

where the second term reflects the presence of a background charge Q = (ei + t2)lg s 
in the corresponding conformal field theory. Indeed one can show that the operators 
Lfc satisfy the Virasoro algebra with a central charge given by Q. Using the formula 
for (f)(x) we arrive at the following expression for the energy-momentum tensor as an 
operator 

g 2 s T(x) = 2e x V W ' {x) - ^ + W\xf + ( Cl + e 2 )W"(x) . 

L — ' X — Za 

i 1 

Then one can show that the Ward identity satisfied by the partition function takes 
simply the form of the Virasoro constraint satisfied by the ground state wave-function 

L n Z = , n > —1 . 

Now we will turn our attention to the brane partition function. We will see that 
also the Ward identity equation for the brane partition function can again be written 
in terms of Virasoro constraints. But instead of the Virasoro ground state condition 
satisfied by the closed partition function, the brane partition function satisfies equations 
analogous to the BPZ equations satisfied by the degenerate states corresponding to 
reducible representations of the Virasoro algebra. In particular, the two types of branes 
of the matrix model exactly correspond to the two types of degenerate states in the 
Virasoro minimal model and Liouville conformal field theory. 
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The comparison with the conformal Ward identity is analogous to the vacuum 
equation discussed above. The only difference is, just as the energy-momentum tensor 
acts on the degenerate insertion as well as other operator insertions, in the matrix 
model energy-momentum tensor we should use the effective potential W(y) perturbed 
by the brane insertion at location x and related to the original potential W(y) as 

W{y) = W{y) - ^ log(z - y) . 

This correction gives an extra term 

— I ^-dy = W\xf - ^W"(x) 
2m J y - x e a 

in the operator 

2m J y — X 

Using the energy-momentum operator, the differential equation (I4.5P can be written 

as 

(-e 2 a d 2 x + g 2 s L_ 2 )Z a (x) = 0. (4.10) 
Notice that, upon a simple rewriting of the above equation as 

(6 2 «9 2 + L_ 2 )Z 1 = 0, (b~ 2 d 2 x + L^ 2 )Z 2 = 0, 6 2 = ^, (4.11) 

the Ward identity equation with a brane insertion takes the form of the BPZ equation 
satisfied by a correlation function in two-dimensional CFT with a degenerate operator 
inserted at location x. 

5 Quantum Riemann Surface and the NS Limit 

In section 14. 1[ we have seen that an interesting limit to consider in the /3-ensemble 
matrix model is when 

ei — > with e 2 = h finite . (5-1) 

In this limit, the ei-brane decouples, as the corresponding coupling constant t\ vanishes. 
However, the coupling constant of the e2-brane remains finite, and the time-dependent 
Schrodinger equation (I4.6p reduces to a time-independent Schrodinger equation 

(-V 2 + W\x) 2 + f(xj\ W 2 {x) = , with [p, x] = H. (5.2) 
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This is the limit considered in [7J in the gauge theory context, and we will refer to it 
as the NS limit. 

This section will be devoted to the study of refined topological strings in the NS 
limit. Notice that in this limit we have g s = — > and we infer that the refined 
topological string theory becomes classical, since g s continues to play the role of the 
closed string coupling constant in the refined theory as we discussed in section [3j This 
constitutes a significant simplification. 

For the more general class of geometries introduced in section [2} we do not have a 
similarly rigorous derivation of the Schrodinger equation satisfied by the brane wave- 
function as for the (subclass of) geometries with matrix model duals. In particular, in 
general we do not have enough input to settle the ordering ambiguities. However it is 
natural to extend the methods of section H] to the /3-deformed Toda matrix model [2J, 
and show that the same result follows. Moreover, note that by taking the rank r of 
the Lie algebra of the Toda matrix model arbitrarily high, we can get an arbitrarily 
high degree polynomial for H(x,p). We leave the details to future work. Here, we 
will just assume that what we derived from matrix model is a general phenomenon. 
Namely, for a B-brane on any Riemann surface S : H(x,p) — arising from a Calabi- 
Yau geometry, as described in section [2j we claim that in the NS limit the e2-brane 
satisfies a time-independent Schrodinger equation with the Riemann surface as the 
Hamiltonian, i.e., 



with ty(x) := ^2{x)- We will test this idea in section [7J for several non-trivial toric 
geometries and observe that the conjecture checks extremely well, despite the absence 
of a matrix model and the presence of ordering ambiguities. A related observation that 
the gauge theory instanton partition function in the presence of the surface operator 
becomes in the NS limit the eigenfunction of the Hamiltonians of the corresponding 
quantum integrable systems has been previously made in [TTl [TT| [91 1^8] . 

In the next subsection, we show how the brane wave-functions transform under 
monodromy in the NS limit. This will be useful to elucidate the observation of [7J 
relating the instanton partition function and the energies of certain quantum integrable 
systems, as we will explain in the following subsection. 





(5.3) 
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5.1 Brane Monodromies and the NS Partition Function 

As discussed in sections 12.21 and I3.4[ the monodromies of the exact wave-functions of 
the branes put constraints on the closed topological string partition function. This 
is true both for the usual as well as the refined topological strings. Combined with 
the Schrodinger equation (15.31) describing the probe e2-brane wave-function, this gives 
important and computable information about the closed string partition function in 
the NS limit. 

Solving the Schrodinger equation for ^O^) and writing 

we conclude that taking the brane around 7# = ^jrrijB 1 induces a change in the 
phase of the closed string partition function by 

Z top (a)-> e^^ dS Z top (a). (5.4) 

On the other hand, as we have seen in section I3T41 on general grounds, taking an e2-brane 
around 7# changes the closed string partition function by 

2 

M~ jB : Ztop(a) ->■ Z top (a + %n) . 

On the other hand, we know that the following is true in the NS limit we have g s — > 
and hence 



Z top (a\ ei, e 2 ) = exp ( ^J" (0) (a; h) 
\9s 



where the dotted terms are suppressed by a factor of g 2 s . On the other hand, the finite 
shift due to the brane monodromy becomes infinitesimal, so Z top (a) changes by 



^lim M lB : Z top (a) exp | - ^ mid ai ^(a; H) j Z top 



a) . (5.5) 



Consistency of the two equations (15 .41) and (I5.5P implies that as we take the 62-brane 
around a B 1 cycle, its phase changes by 

/ OS = d ai J r{0) (a;h). (5.6) 
Jb 1 

Note that this is accompanied by 

/ dS = ai {h), (5.7) 
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by definition of what we mean by the Aj-cycle. Solving for the brane wave-function 
explicitly from (I5.3P given the classical curve equation X : H(x,p) = (involving some 
arbitrary coordinates on the Calabi-Yau moduli space) and computing the periods ( 15.61) 
and (15.71) from brane monodromies around the corresponding cycles hence defines a 
'quantum' special geometry on the Calabi-Yau moduli space. This is an ^-dependent 
generalization of the usual special geometry. In particular, the aj(h) periods obtained 
in this way are good, 'flat' coordinates on the moduli space. 

5.2 An Explanation of the Results of Nekrasov-Shatashvili 

From the point of view of four-dimensional gauge theory, the NS limit (15. ip is clearly 
an interesting one. In this limit, two of the four dimensions where the Poincare invari- 
ance is broken become effectively compactified and we are left with a system that is 
effectively two-dimensional. Indeed, in [7] Nekrasov and Shatashvili observed that the 
gauge theory in the f2-background in this limit is connected to various quantum inte- 
grable systems in a highly non-trivial way, a feature shared by various two-dimensional 
supersymmetric gauge systems [49], [50], [5T] [52] [53] [26] . The Bethe ansatz equation for 
the corresponding integrable system turns out to coincide, according to [7J, with the 
critical points of the free energy in this limit. More precisely it is given by 



The function W is related to the refined gauge theory partition function Z/v =2 (a; ei, £2) 
at arbitrary via 



Furthermore, the corresponding eigenvalues of the g commuting Hamiltonians can be 
expressed in terms of functions of the solutions aj to the above equation. The question 
to address is why studying the partition functions of a 4d gauge theory system and 
taking such a limit should have anything to do with answering these questions for an 
integrable model. Here we would like to explain this fact. The basic point will be the 
following: We can study an M = 2 gauge theory system in 4d in two different ways. 
The first one uses a target space description while the second one uses world-sheet 




(5.8) 



These equations determine aj given a g-tuple of integers n 1 by 



— = 2mn . 



W = lim ei log Zj^ =2 . 
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techniques. The fact that these two approaches have to render the same results will be 
the key to an explanation of the observation in [TJ. A related consideration has also 
appeared in [12]. See also [53] . 

In particular, the NS computation is the target space viewpoint, and the integrable 
system emerges from the world-sheet viewpoint. The equality of Zj^ =2 and the refined 
topological string partition function Z top implies 

W(a,h) = ^°\a; ft), 

where J-^^a; ft) is the genus zero refined topological string amplitude. 

Consider an M = 2 system given by a Seiberg-Witten curve E : H(x,p) = 0. We 
assume, as in the rest of this paper, that the corresponding SW differential is A = p dx. 
This theory can be engineered in type IIB strings by considering a local CY given by 
a hypersurface with the equation 02 .ip . or in M-theory in terms of M5 brane wrapping 
the curve E in the complex two dimensional space (x,p) C C 2 . Let us assume this 
curve has genus g. Then we know that there are g deformations of this curve given 
by the g Coulomb parameters aj=i,.., s . Let us further assume, as is the case in all the 
known examples, that the a\ dependence of H can be written as 

9 

H(a;x,p) = H (x,p) + f I (x,p)E I (a) 

i=i 

for some parameters Ej that are functions of the flat coordinates a/. The NS pre- 
scription turns out to translate in this language to the statement that aj are fixed 
by the 'critical points' of W and these specify the eigenvalues Ej of the g commuting 
Hamiltonians. 

The key idea to connect the integrable system to the refined topological string 
considerations is a trick known as the 'separation of variables' (SOV) [53]. Think 
about C 2 as a complexified classical phase space. Consider g points on it given by 
(pi,xi), as / runs from 1, ...,g. We can choose the moduli of the SW curve to pass 
through all these points. In particular this fixes all the Ej in the above equation in 
terms of these g points: 

Ei = h I (p 1 ,x 1 ;p 2 ,x 2 , •••) , 
for some functions hj that are readily computed from requiring that 

H(pj,xj) = , J = !,...#. 
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It is not difficult to see that the defined hj are classically, as well as quantum mechan- 
ically, commuting and thus define an integrable structure (i.e., we have g commuting 
Hamiltonians hi for a 2g dimensional phase space (pi,xi)) [56]. Furthermore, it can 
be shown that in the quantum theory the eigenstates of hj can be written using the 
separation of variables in terms of single-particle wave-functions: 



where each ^>j(xi) is annihilated by H(x,p). Namely, they separately satisfy the equa- 
tion (15.31) . In the context of refined topological strings, we identify the corresponding 
wave-functions *f*(x) as the amplitudes for the branes in the refined open topological 
strings, in the NS limit. We are almost finished except for the fact that the brane 
wave-functions are holomorphic functions defined on the complex x-plane and not just 
on a real space. Hence, the wave-functions could have monodromies and the question 
of their single-valuedness arises. To study this issue, let us assume that we have chosen 
a real subspace as well as a choice of the real coefficients hi in the SW geometry, such 
that H = leads classically to g physically allowed intervals in the x space. Let us call 
these intervals b 1 , for 7 = 1, ...,g. Note that the end points of the classical intervals 
correspond to p = and in the hyperelliptic case, and correspond to branch points of 
the SW curve over the x-plane. Let us assume that the exact wave-functions have n 1 
zeros in the interval b 1 . Consider cycles B 1 in the complex plane encircling b 1 . These 
would correspond to a basis of B-cycles of the SW curve. The fact that \& has n 1 zeroes 
on the interval b 1 implies that in the analytically continued picture, surrounding the 
B 1 cycle picks up a phase exp(27rm 7 ). This would also guarantee that the restriction 
of to the real slice is single-valued and well-defined. On the other hand, we have 
seen in the last section in ( 15. 5 p that taking the brane around the B 1 cycles changes its 
phase by the derivative of the closed string partition function di\V(a): 



Putting the two together, we conclude that 



thus explaining the results of NS. This relation between the phase of the wave-function 
and the integers n 1 is known as the 'exact Bohr-Sommerfeld' quantization rule. It 
generalizes the usual rule to a form which is true to all orders in h. In this context, 
the relevance of the NS limit was also pointed out in [57J EHl E] • 



<ff(x l ,...,x g ) = JJ 



i 
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The scaling n 1 ~ (3N 1 

The quantization of n 1 for the solutions of the Bethe-ansatz equations is reminiscent 
of the filling fractions in the context of matrix models. Let us consider the cases where 
the topological string theory is dual to a matrix model. Then consider a closely related 
matrix model with a number of critical points, for which the above physical intervals 
b 1 become the cut around the critical points in the eigenvalue plane C at large N. By 
considering this other matrix model we are effectively exchanging what we call the A- 
and the B-cycle. Denote by N 1 the filling fractions for each of the critical points. As 
noted before, the NS limit corresponds to taking e± — > 0, N 1 — )■ oo while keeping the 
't Hooft coupling ft 1 = eiN 1 finite. And the period around the cycle B 1 encircling the 
cut is given by 



Jbi 

On the other hand, in the WKB approximation the Bethe-ansatz equation reduces to 



Of course this relation is only valid to leading order in h = 
6 Genus Zero Examples 

In the context of matrix models we have shown that if we consider the Gaussian 
matrix model the brane wave-functions satisfy the exact equations H(x,p)^/ a (x) = 
with [p,x] = e a , where 



In particular the time-dependent Schrodinger equation (14.61) has trivial time depen- 
dence that has been absorbed in the shift of the background, and therefore we obtain 
the same equation in the NS limit as in the general case without taking the NS limit. 
For multi-cut matrix models with have higher degree potentials, this is no longer true. 
In these cases the genus of the spectral curve is greater than zero, and the multi-time 
dependence of the Schrodinger equation captures the moduli of the higher genus curve, 
which are absent for the genus zero case. It is thus natural to conjecture that for 





Thus in this limit we have 



n 7 e 2 = N I e 1 n 1 = N 1 (3 . 



H(x,p) = —p 2 + x 2 — jJL . 
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all the cases where S : H(x,p) = is a genus zero curve, as in the Gaussian matrix 
model example, the time dependence of the Schrodinger equation is absent. As already 
noted, the exact equation satisfied by the branes for any e±, e 2 is the time independent 
Schrodinger equation. 

In more details, on a genus zero Riemann surface, in general we have a set of 
compact A-cycles. All the B-cycles are by necessity non-compact, though we can find 
a basis of them satisfying Aj fl B J = 5j. The monodromies around the A-cycles define 
the flat coordinates 



In particular, the partition function of the brane depends only on one e a , associated 
to the corresponding brane 



Solving the Schrodinger equation, we can find not only the brane partition functions, 
but also the closed string partition functions. In general, on a genus zero Riemann 
surface we have a set of compact A-cycles. All the B-cycles are by necessity non- 
compact, though we can find a basis of them satisfying AjC\B J = 5j . The monodromies 
around the A-cycles define the flat coordinates 



Consider now taking the brane brane around a S-cycle. On the one hand, from 
general considerations 



where J-"(a; ei, e 2 ) = g 2 log Z top is the exact free energy of the refined topological string. 
On the other hand, since the solution for the brane is independent of e 2 , and taking 
the NS limit e 2 — > 0, we find that the same amplitude equals djW: 




^ a (x) = exp(S(x,e a )/e a ) . 





dS/e 1 = T{ ai + gl^/g] - F( ai - gl/le^/ g 2 s 



S 1 




Hence, we derive that 



F{ ai + gl/Ze^/gl - F{aj - gl/^/g 2 = d T W , 
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and both equal to the change of phase of the brane, after being transported around 
the B-cycle. The fact that the two expressions agree is special for genus zero Riemann 
surfaces. In general, monodromies of the branes give rise to difference equations for 
the partition functions. In the NS limit, the difference equation becomes a differential 
equation, for W. 

In this section we present two classes of examples along these lines. The first one 
involves the Gaussian matrix model and the second class involves toric geometries 
without compact 4-cycles. 

6.1 The Gaussian 

The simplest matrix model to which our derivation of the Schrodinger equation of 
section 0] applies is the Gaussian 



which has one A-cycle, corresponding to the circle at real values of x and p, and one 
non-compact B-cycle. The /3-ensemble matrix model partition function can be easily 
evaluated explicitly to find the free energy (under the shift fi — > [i + (ei + 62)/2) 



One infers that this is exactly the partition function of the c = 1 string at radius 
R = [2], which has been computed in [SS]. Note that the NS limit in this case 
corresponds to the R — > oo limit at fixed [i. In particular, the partition function per 
unit volume corresponds precisely to the NS definition of the partition function in this 
limit: 




i<j 



The classical spectral curve is the deformed conifold 



H (x, p) = —p 2 + x 2 — /i = , 



(6.1) 




W(/i) = lim ei log Z , 



with 




Note that from this it follows immediately that 
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in agreement with our claim. We want now to check that we can indeed recover 
in the NS limit the coefficients (j6.2p via the monodromies of the brane solving the 
Schrodinger equation we derived in sections H] and El For this, we will need the power 
series expansion of W: 



where 




Adding an ei-brane, the Riemann surface becomes quantum: 

This is the (inverted) quantum harmonic oscillator. The monodromies of the wave- 
function ty(x) around the B-cycle have been in fact calculated to all orders in h already 
some time ago in the context of Stokes phenomena in quantum mechanics, see for 
instance [60J. (We let h = e 2 for the rest of this section.) The result can be matched 
with (16. 2p . Alternatively, they can also be computed using the matrix model techniques 
of solving the loop equations. This method has been used to solve the /3-ensemble 
Gaussian matrix model in [M] . 

In order to illustrate the technique we will use later in the context of the more 
complicated cubic potential in section [7J let us re-derive this result in a simpler, but 
non-exact fashion. In detail, we can solve the Schrodinger equation for the brane wave- 
function ty(x) via the usual WKB Ansatz known from elementary quantum mechanics, 
i.e., 

V(x) = exp (^S(x)j , (6.3) 

with 

00 

S(x) = S (x) + Y,Sn(x)h n . (6.4) 

n=l 

Plugging this Ansatz into the Hamiltonian and expanding in h yields for the first few 
orders 

oc o p l^i Ac (d x \ogp) 2 -2d 2 x \ogp 

d x S = p , d x S ± = --d x logp , d x S 2 = , . . . . (6.5) 

2 Sp 

Note that the WKB approximation determines directly not S(x), but d x S. Near any 
point on the Riemann surface, we can integrate this to obtain S(x) locally, however 
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the global solutions have monodromies. In particular, to obtain from the semi-classical 
wave-function ty(x) the (quantum corrected) periods lie = Yl^o^c^ ft 1 , we have to 
integrate dS along the classical cycles of the geometry. 

For specific examples, or better potentials V(x), things considerably simplify, since 
we can turn the contour integration for the higher order corrections into a differentia- 
tion. For example, this idea has been put forward in the gauge theory context in |57j . 
In detail, we can argue for differential operators T>^ such that 

n (n) = V {n) n (0) _ (6 6) 

Defining 



n=l 



the full quantum corrected periods are simply given by n c = vn K c 



( o) .Q 



We can partial 
nient expressions 



y integrate the contour integrals of ( 16.51) to obtain the more conve- 



S' (x) = - E) 

V " [X) 



A%{V{x) - Ef/ 2 ' 

o, (r] _ WW 2 V""(x) 

4W " 1536(y(x) - E) 7 / 2 768(F(x) - Ef/ 2 



(6.7) 



where we used the parametrization p(x) = \fW' [x) 2 + f(x) =: y/V{x) - E. For the 
Gaussian, we can identify V(x) = x 2 and E = fi and the operators T>^ are particularly 
simple because V"{x) = 2 and d> 2 V(x) = 0. Namely, one has 

V^ = --d 2 u , V^ = +—d A u , V^ = ^— «9<\ .... (6.8) 

24 5760 967680 K ' 

We conclude that the period IT^ does not receive any quantum corrections, while we 
obtain for the B-period Hb = D 

4 Note that we know from the M-theory picture described in section f3. II that the refined partition 
function in the NS limit is odd in H, hence 2?( n ) with n odd vanishes. In particular, since S' n is always a 
total derivative for n odd, this implies the integral of total derivatives vanish and there is no subtleties 
involving singularities. 

5 For simplicity, we absorbed a normalization factor of \/2. 
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U B (a) = n ( R 0) (/i) - — u~ l h 2 + —u~ 3 h A - -^—u~ 5 h 6 + . . . , (6.9) 
bkh>j b vav 24 ^ 2880 P 40320^ ' V ; 

or, after integrating over /x, 

fi 2 7 31 

= T (0) (u) - — log// — /i" 2 fr 4 + — — /i~ 4 n 6 + . . . , (6.10) 

vw vw 24 5P 5760 P 161280^ ' V ; 

in agreement with the expectation (16. 2p . 
6.2 Toric Geometries without 4-cycles 

Consider the A-model topological string on a toric Calabi-Yau X without compact 
4-cycles. This is mirror to the B-model on the Calabi-Yau 

uv + H(e x , e p ) = 0. 

The corresponding Riemann surface has genus zero and is of the form 

H(e x ,e p ) = P n (e x ) + e p+mx P k (e x ) = , (6.11) 

where P n ,k{e x ) are polynomials of degree n and k in variable e x , and m is an arbitrary 
integer. This Calabi-Yau has n + k — 1 moduli which enter as coefficients of the 
polynomials. The simplest example is the conifold with n — h — 1 and m arbitrary. 
As explained in [T3| WI\ , there is a particularly natural way of writing the Calabi-Yau 
geometry in terms of open and closed string flat coordinates that are typically given 
by certain periods on the mirror geometry. By a change of variables we can rewrite 
(I6TTT|) as 

71 k 

H(e x ,e p ) = - Q ai e x ) - e p+mx ]^[(1 - Q^e x ) . (6.12) 

1=0 J=0 

Depending on the chamber, it is natural to set one of Q ai or Qp 3 to 1, which we can do 
by a leftover degree of freedom to shift x and p by constant values. This gives indeed 
altogether k + n — 1 moduli. The Riemann surface can be viewed as a copy of a u 
cylinder, where the one-form 

d(p = pdx 

has singularities at x — t ai , x — t$j. Up to a gauge degree of freedom discussed above, 
x, t ai = —\ogQ ai and tp s = — logQ/3 7 are fiat coordinates on the open and closed 
moduli space. The canonically conjugate variables are 

[p, x] = ei 
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In particular, the brane satisfies the exact Schrodinger equation 

H(e*, e p )V(x) = , (6.13) 

which is a difference equation, 



A' 



- Q aj e*gr 1/2 )*(x) - JJ(1 - Qp^q^Wx + H) = , 

7=0 J=0 

where gi = e _ei . Note that we have redefined the Kahler parameters £ a and tp with a 
shift. This has exact solutions in terms of quantum dilogarithm functions 

where 

oo 

L(e x ) = - q[- l/2 e*) . 

Note that L(e x ), given by the quantum dilogarithm, satisfies 

L{e x+ei ) = (1 - q; 1/2 e x )L{e x ) . 
This includes, for example, the conifold, where 

* (r) - L{6&Q) 

^ conifold\ u 'j T i x.\ 

L[e x ) 

solves (I6.13P with 

H comfold (x,p) = (1 - eV l7 V " (1 - e x Qq; 1/2 ) . 

Note that in general there are ordering ambiguities in defining the Hamiltonians, but 
they can all be absorbed into the shifts of the open and closed moduli, i.e., what we 
mean by x, t. As a consequence, in this particular case, there is really no physics 
behind the quantum shifts, and we can simply choose them in the way that is the most 
convenient. Any two such choices differ by re-parametrization of the moduli space, and 
hence are physically equivalent. Finally, the above expressions all assume we are in the 
regime of the moduli space where Qe x < 1. If that fails to be the case, for example 
due to moving the brane around the Riemann surface, it is better to use the analytic 
continuation 

L(e x ) = e -(* 2 +i™)/^L(e- x ) 
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to rewrite the wave-function. 

Consider now moving the branes around. First, it is easy to show, in these coor- 
dinates the variables t and x coordinates we defined via ( 16. 12ft are already the open 
and closed A-periods. Secondly, since there are no compact B-periods, the branes are 
essentially single-valued on the Riemann surface (the only shifts come from the classi- 
cal piece of the wave-function). As a consequence, computing the B-periods amounts 
to evaluating \1/ at the poles of the Riemann surface. The good B-periods are obvious 
either from the A-model - they are mirror to non-compact 4-cycles, or equivalently, by 
asking that the only infinities in evaluating the periods come from the classical pieces 
of the wave-function, corresponding to the fact that at infinity, all instanton correc- 
tions are suppressed. Consider, for example, the B-period corresponding to bringing a 
brane in from infinity to a singularity with x — t ai , and sending the brane back out to 
infinity at the singularity with x — > oo. Correspondingly, we find 

exp (I f ^ - *<*«" 



ei J B i J *(oo) 

up to the contribution of the classical pieces x 2 je\ to the answer - these are ambiguous 
and reflect the non-compactness of the geometry. On the other hand, this non-compact 
B-cycle intersects all the 2-cycles whose areas are of the form 



J A j 



by 

B 1 n A = n 1 . 
This implies that the free energy jumps by 



e i J b 1 J Z(t ai — y 
or 

tf(t aj ) _ Z(t ai + f ) 



tt(oo) Z{t ai - f ) ' 
To compute ^>(t ai ) we need to analytically continue all the terms in the solution above 
corresponding to which Qe x becomes greater than one at x = t ai . There are similar 
equations with t a replaced by tp. We can view this as n + k equations for the partition 
function Z(t a ,tp) depending on n + k parameters, which we can solve, with the bound- 
ary condition, that in the limit where all the t's we recover the partition function of 
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n + k disconnected copies of C 3 . It is easy to see that the solution to these equations 
is the refined partition function of the closed topological string on this Calabi-Yau 

Z(t t ) = M(l) n+k ^o</<J , <n M (Q<*io<i>) Uo<j<j'<k M (QpjP., i ) 

\ Q / ' Pj) \ ) i—\n,k 



where, e.g, Q ai p,, is understood to be QaJQpj if Re(t a/ ) > Re(t^ 7 ) and Qpj/Q ai 
otherwise, and similarly for the rest, and 

oo 

m(q)= n (i-vr^t^Qr 1 



is the refined MacMahon function, where we have defined q 2 = e e2 . In the above we 
have repeatedl 
we get simply 



have repeatedly used relations such as L(Q) = m ^q -1/2^ ■ For example, for the conifold 



z (t) mi] 

n conifold\ b ) 



M{Q) ■ 

The fact that this is also consistent with the monodromies being computed by W 
can be seen as follows. Note that we can write the free energy of the topological string, 
Tjgl = log(Z) as 

n,k 

t IqI= Y 7(w*,/)+ Y ittwj') - Y 7 (Wj)> 

0<I<I'<n 0<J<J'<k I,J=0 

where 

7 (t;e 1 ,e 2 ) = log M(Q) . (6.14) 

In particular, 

l(t, ei, e 2 ) = V , 

with [n] ea := (qa 2 — q a n ^ 2 ), and Q = e~*. On the other hand, the reduced partition 
function W = lim^^o J 7 1 e i is built out of the function 

Qn 

lNs(t,e 2 ) = 



i n 2 [n] £2 

n,k 

>V= Y lNs(ta iai ,)+ Y lNs{tpjPj,) ~ Y 7ATs(Wj) ■ 



0</</'<n 0<J<J'<k I,J=0 
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The fact that 




exp(7 dS/e 2 ) =ex P (d ta W), 



\Jbi 



in agreement with our claim. 

7 Genus One Examples 

Let us now consider some more non-trivial geometries and confirm that also for them 
the NS limit of the refined topological string partition function can be recovered from 
a brane probing the quantum geometry, as expected from the discussion of section |5j 



The /9-ensemble matrix model partition function can be calculated perturbatively, as 
for example done in [63]. The most leading ^-correction to the free energy of the 
cubic /3-ensemble matrix model has also been computed in [39] using the loop equation 
techniques. Here we will use the techniques discussed in section and obtain more 
non-trivial higher order results. 

As a check of the statements of section 15. 1[ we now want to compute the first few 
terms in the ^.-expansion in the refined topological string free energy for the cubic 
by studying the brane wave-function satisfying the Schrodinger equation of section [51 
which can be obtained by simply canonically quantizing the spectral curve. For that, 
note that in the case of the cubic, the 1-form of the matrix model spectral curve (I3.10P 
can be expressed in terms of the branch-points (roots) 2?i,...,4 as 



7.1 The Cubic 



Let us consider the matrix model with cubic potential 




A = dx g\f( 




(7.1) 



Following [35], we define new variables 



Q = o 
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In terms of these variables, explicit expressions for the periods , Iljg , ngj and 
can be easily found via expansion for small z\ and z 2 of the corresponding contour 
integrals of A (see appendix B of [35]). As usual, we will set Q = —5 and X = 
y/5*-2( Zl + z 2 ). 

On the other hand, it is known that the B-model refined 1-loop amplitude on this 
geometry is given by 

1 3 

7^( Sl , s 2 ; f3) = - logdet g + ]T kj{(3) log(Aj) , (7.3) 

7=1 

where denote the flat coordinates given by the A-periods, are other coordinates 
given in (17. 2p . gij = d Si Zj, A/ are the three discriminant loci (c/., [64] ) 



Ai = Zl z 2 , A 2 = ^5 2 - 2( Zl + z 2 ) , A 3 = (5 2 - 3 Zl ) 2 - Qz 2 5 2 + 9z 2 2 + U Zl z 2 , 

where kj((3) are /3-dependent coefficients. In particular, we expect that in the B-model 
all the information about the refinement is encoded in the coefficients Ki(/3), hence is 
of a purely holomorphic nature. 

It remains to fix Ki{f3) by the behavior of the topological strings amplitude near 
the singular points. Since Ai is a conifold like singularity, we expect from [1] that 

«i03) = ~{P + r 1 ) ■ (7.4) 

The other two k 2 ^{f3) can be fixed by considering the special moduli S\ = —s 2 = s(8) 
that is related to the SU(2) Seiberg-Witten theory [65]. By matching the refined 1-loop 
amplitude at these special moduli to the £77(2) partition function in the ^-background, 
we conclude 

«2G9) = ~, « 3 (/3) = ^(6 - - P- 1 ) ■ (7.5) 

Comparing with the explicit perturbative matrix model expansion shows that (17. 3p 
with (I7.4p and (17. 5p indeed reproduces the refinement of the cubic results, and hence 
we have verified that the refined B-model is indeed dual to the /3-ensemble, at least at 
the 1-loop level. 

For later comparison let us give the following explicit expansions of J-^^si, s 2 ; /?), 
obtained using ( 17.30 and (17.40 and (17.51) 

lim d 2 F {1 \s, -s; 0) = —s~ 2 + 153 + ^^-s + 1217160 s 2 + . . . , (7.6) 

/3^o s v ' 12 3 
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lorn p d 2 si s 2 ; (3) = -sf + - + — Sl - — s 2 — Sl s 2 + . . . , (7.7) 

where we have set for simplicity 6 — g — 1. @ 
Seiberg-Witten slice 

Let us consider first the simplified case with s = s\ — — s 2 (the Seiberg-Witten slice). 



This restriction of parameters has been discussed in detail in [65] . The one- form takes 
the simpler form 

A = dx \/ {x 2 — a 2 ){x 2 — b 2 ) , 

and the periods can be easily expressed through complete elliptic integrals of first and 
second kind, i.e., 

< = f A = A ((a 2 + b 2 )E(ki) - 2b 2 K(h)) , 
2m J a bir \ J 



ng } = 2 J" A = ^ ((a 2 + 6 2 )£(fc 2 ) - (a 2 - b 2 )K(k 2 



(7i 



with /c 2 = 1 — \ and = The relation between the parameters a, b and the original 
parameters 5, A reads 



a = -^VS 2 + 8A 2 , 6 = -^V5 2 -8A 2 . 

Inverting LT^ (A), plugging it into 11^ (A) and taking the derivative yields the 
gauge kinetic function of the cubic matrix model under specialization si = — s 2 = s 

[SSJ- 
Let us now apply the approach of section [5l using the technique introduced in 
subsection 16.11 It is convenient to absorb the energy E into the potential V(x) such 
that V(x) = (x 2 — a 2 )(x 2 — b 2 ) and E = 0. Hence, at order h 2 we have to integrate 
(see (J62D) 

dx YM (7 9) 

48((x 2 -a 2 ){x 2 -b 2 )fl 2 ' V ' ' 



1 



6 Notice that, as explained in the previous sections, from the quantum mechanics point of view it 
is more natural to consider a double expansion in g 2 = — e\C2 and H = ei + 62, while from the point of 
view of the holomorphic anomaly equation it is more natural to expand in g s and treat ft = — ei/e^ as 
a parameter. Also notice that the parameter j3 always enters in the symmetric combination j5 + 
in the latter formalism. In order to compare the results obtained from these two viewpoints, simply 
notice that g 2 ((3 + /3 _1 ) = h 2 + 2g 2 . So the coefficients lim^o (3 J 7 ' 1 ) should really be compared with 
the g®h 2 term in the g s , h expansion when taking the NS limit. More generally, the term lim ( g_ i> o P 9 
contributes to the term g®h 2g in the g s ,?i-expansion. 



16 



It is convenient to define a differential operator T)^ acting on 

1 1 
dx — — = K(k\ 



b 1 1 

/ ax — = -fM^'i, 

2niJ a ^( X 2 - a 2 ) (a; 2 - b 2 ) 2ira 

and 

b 1 2 

dx — = = -K(ko) , 

6 ^(x 2 -a 2 )(x 2 -b 2 ) a 

to obtain ( 17.9ft . It is easy to see that 

^ (2) = ~~pTTT~2 — Tz\ W ~ 5 ^ ~ b (b 2 - 5a 2 )d a ) , 
bab{a z — b z ) v 7 

does the job. Hence, (I7.9P evaluates to 

S {2) = - J^a^W) (( ° 4 " 10a%2 + bA)E{h) + Ab2{a2 + h2)K{kl)) ' 
n(2) = - o h 2 ( \ ixs - 10a%2 + b ') E ( k 2) ~ (« 4 - 6aV + 5b A )K(k 2 )) . 

Proceeding for s := 11^ + h 2 U^ and IL3 = + ft 2 n^ as sketched for the pure 
classical part above, yields the first order quantum correction to the gauge kinetic 
function. The result precisely matches the expectation (17.61) . 

General case 

Let us check that the agreement also holds for the general cases where si and s 2 are 
independent. The starting point is as in the previous section the one- form expressed 
through the roots given in (17. ip . We have to integrate 

dx K!W (710) 

48((x — x\)(x — x 2 )(x — x 3 )(x — X4)) 3 / 2 

where now V(x) = (x — X\)(x — x 2 )(x — x 3 )(x — £4). As above, this can be achieved via 
defining an operator acting on the classical periods. In order to construct such an 
operator, note that it is sufficient to find operators Vj with / = 0, 1, 2 satisfying 

VjV(x) = x 1 , 

since V"(x) is a polynomial of degree 2. The operator T>^ is then a combination of 
the T>j. Such a basis of operators for the given potential can be inferred to be 

-(x 2 - x 3 )d Xl + (xi - x 3 )d X2 - (xi - x 2 )d X3 



Oi - x 2 ) (x 3 - x 4 ) V ( x i ~ x 3 )(x 2 - x 3 ) 

(x 2 - x 4 )d Xl - Qi - x 4 )d X2 + (xi - x 2 )d X4 
(xi - x 4 )(x 2 - £4) 

17 



„ 1 / (x 2 - £ 4 )(3e 2 - c^J - (xi - x 2 )(c) :E2 - d X4 ) 

L>i — X 4 - 



+ X 



(x 4 - x 3 ) V (xi - x 2 )(x 2 - x 4 )(xi - x 4 ) 



(xi - x 2 )(xi - x 3 )(x 2 - x 3 ) 

1 4 x 

22 n ~ x K)x]d XI 



nJ<jO/ - a; J 



The combined operator reads 

P( 2 ) = - (c 2 D 2 + dl?! + C P ) X>0 , 

with 

c = 2(x 2 x 3 + (x 2 + x 3 )x 4 + x 4 (x 2 + x 3 + x 4 )) , 
ci = -6(xi + x 2 + x 3 + x 4 ) , 
c 2 = 12 . 

Applying T)^ to the classical periods yields the first quantum corrections. Inverting 
Si := (1 + /i 2 l}( 2 ))n^ yields the quantum correction to the mirror map at order h 2 , 
i.e., 

s 2 )| fl 2 = 2 + 62si - 78s 2 - 4808sis 2 + . . . , 
22(^1, s 2 )| fi 2 = 2 + 78si - 62s 2 - 4808sis 2 + . . . . 

Plugging these into n 4 = (1 + h 2 T>^)U^ and taking the s\ derivative yields the 
expected expansion (17.71) . 

7.2 Local P 1 x P 1 

Let us consider non-trivial toric Calabi-Yau geometries of genus 1. For example, the 
mirror curve for local P 1 x P 1 can be parameterized as 

H(p,x) = -1 + e x + e p + Zl e~ x + z 2 e~ p = . (7.11) 

The skeleton of the corresponding Riemmanian surface is illustrated in figured] together 
with the integration contours for the classical periods [66] . 

Inspired by the matrix model results of section we consider the brane wave- 
fucntion satisfying the difference equation 

(-1 + e x + zie~ x )m(x) + *(x + h) + z 2 ^(x -h)=0, (7.12) 
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obtained by canonically quantizing the curve (17. lip . In the limit x —> oo, we can view 
(I7.12p approximately as a difference equation with constant coefficients, whose solution 
takes the form 

V(x)=ir 1 (x)e S £ x + Tr 2 (x)e ! % x , (7.13) 

with some functions 717(3;) periodic under x — > x + H and some constants a/. In 
particular, for a x > a 2 one can show that [6] 



e- = li m + . (7.14) 



However, comparing (IT. 13[) with the general form of the brane wave-function near 
x — > 00 which can be inferred from [62], one infers that actually 

ai = - (U Al - log zi) , 1 = 1,2, 

where 11^ denotes the fully quantum corrected A-periods. Explicit calculation of a± 
gives 

ai = - Oi +Z2)- ( 4 + ^ + g ) Zl Z2 ~~ \^ + ^2) 

1 O 6 \ , o O x 10 



16 + + q 2 + - + Qq) {z x z\ + z\z 2 ) - —{z\ + z\) + . . . 
q q J o 

as the first few orders in small 2^2 expansion, where q := e h , and we have kept only the 
finite terms. Expansion in h yields 11^ = ^2^L h n with 11^ the classical period 
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and the first few quantum corrections 



n (1) - 

= -2z 1 z 2 - 20(zfz 2 + Zi4) ~ 420^2 + . . . . 
(3) (7-15) 

nf = o . 



Inverting gives the mirror maps including quantum corrections thereof. We obtain 

zi(Qi,Qa)Ui = zi(Qi,Q 2 )\n : i = 0, ^ 
zi(Qi,Q 2 )\v = -2QjQ 2 - AQlQl + . . . . 

with Qi = e~ li and tj are the Kahler parameters of the geometry. One can also solve 
the difference equation order by order in h using a WKB Ansatz. For the leading 
(classical) piece we obtain 



S' Q = p (x) = log 



-1 + e x + z x e~ x ± e - x y/(-e x + e 2x + z x ) 2 - 4e 2x z 2 



where the two different signs correspond to the two solutions *f>±(x). The higher order 
terms S' n can be easily obtained, though they are too lengthy to be displayed here. It 
is convenient to use the coordinate x = e x instead of x. Let us fix a branch, say ^+(x) 
and let us expand S' n for small Zj. Effectively, this means that we restrict to a local 
C 3 patch of the geometry, as also shown in figure [TJ We observe that in this limit the 
expansion of S' n has poles at x = and x — — 1. The poles correspond to the two 
'internal' punctures of the local C 3 patch. Taking for instance the residue at x = 
yields the instanton part of the classical A-period, i.e., 



^(llig-log*j) = j> ^S' = -(z 1 + z 2 )-Qz 1 z 2 -^(z 2 + z 2 2 )-30(z 2 1 z 2 + z 1 z 2 2 ) + ... . 

(7.17) 

Integrating the higher order SL we can indeed reproduce our previous results f)7.15p . 
as explicit calculation shows. L 7 ! 



7 Dcpcnding on the actual parameterization of the curve, one might have to take the correct com- 
bination of 'small' periods ai, i.e., IT^ = (j£ ai — § a ^j i§ S' n [62], to ensure that the odd sector in h 
drops out. The contour integral around cti can be obtained by performing a SX(2,Z) transformation 
of the curve such that in the limit z\ , zi <C 1 we end up in the corresponding local patch. 
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Let us now consider the B-period. We observe that we can recover the instanton 
part of the classical B-period in the local patch via the integral 

A di 1 2 11 2 77 2 47 2 1377 2 2 ^ 

On — — Zo — — Z-\ — Zn — DZtZo — Zl Zr, — Z-i Zo — Z-i Zn -\- . . . 

s x 2 1 4 2 2 2 2 1 4 1 2 ^ 

up to a constant, and where we have taken 5 — > 0, A — > oo. Note that this integral 
diverges as 5 — > and we kept only the finite terms. Hence, 

f f s »=K <) "i <,iog(2i) "J ios(2i)2 )- 

Similarly, we obtain for the order H 2 

n (2) _ 1 8 „ _ *2 62 2„ 107 2 E fi4 -2„2, 

ll Bl — ~q Z 2 - ^ 1 2 - 2 Z 2 §~ 1 2 3~ 2 _ ^ mZ l Z 2 + • • • • 

Plugging the mirror map (17. 16ft into 11^ = + Ll^/l 2 yields at order ft 2 

n Bl (Qi, Q2)W = c - % - ^ig 2 - f - ^g 2 g 2 - ^g 2 . + . . . , 

o 3 o 2 d 

where c denotes some constant. 

Let us compare this result to what we expect. In the B-model, it is easy to see that 
the refined 1-loop amplitude 

= -log(det^) __l og (A) — log(ziz 2 ) , 

with gjj = \im h ^ dQjZj\ and A = 1 — 8(zi + z 2 ) + 16(zi — z 2 ) 2 reproduces the refined 
vertex results of [3J. Taking the NS limit, we conclude 

U Bl (Q 1 ,Q 2 )\ h 2^Q 2 d Q2 lim/3jrW, 

/3-s-O 

up to the constant part of ^p. 

Let us also check the order H 4 . For simplicity, we will focus on the special moduli 
slice of the moduli space with z = Z\ — z 2 . The mirror map at this order comes out to 
be 

O 3 14 
z(Q)\n* = -^r-^Q 4 + .... 

O 6 



8 The reason why we do not get the constant part in our calculation is a technical subtlety due to 
our method of B-cycle integration rather than a conceptual issue. 
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Figure 2: The skeleton of the B-model geometry of local P 2 . A-period contours are drawn in 
red while B-period contours in blue. Case (i) shows the effective geometry of the parametriza- 
tion (j7d8j) under the limit z < 1, while (ii) of ([772T]) . 



We obtain for the B-period 



Hence, 



„ia\ 33 o 1280 n 

ng> = -2z - —z 2 - —z 3 + 



Ufi W, " 360 6(P 20 ^ + ---- 



Comparison with the refined vertex results of [3] shows that this indeed matches the 
f3 — > limit to the corresponding order under the identification Q := Qi = Q 2 . 

7.3 Local P 2 

Let us consider local P 2 as our final example. We parameterize the mirror curve as 

H (p, x) = - 1 + e x + e p + ze h/2 e~ x ~ p = . (7. 18) 

The corresponding geometry is sketched in figure [2j A canonical quantization yields 
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the difference equation 

(-1 + e x )^(x) + (a: + K) + 2e- x e ri/2 ^(x - A) = . 

As in the previous example, we solve for a>\ defined in (17.141) to extract the quantum 
A-period given by ai = ~ (11^ — log(z)). We obtain 

1 + 9 ( 1 7 7 ,\ , 
ai = - 6 + - + — + -q + q 2 



3 + 9q + 36g 2 + 88g 3 + 144g 4 + 144g 5 + 88g 6 + 36q 7 + 9q 8 + 3q 9 . 



z° + ... . 



3q 9 / 2 

Hence, via expansion in h we infer for 11^ = ^^L 11^ h n the first few corrections 



nW = o, 

U f = J-z - ^ - 630z 3 - 17325^ 4 H 
4 2 



n?=0, (7.19) 

f 4 > 39 „ 2961 , d 
4 4) = - u - - —z 3 - 19635/ + . . . . 



Inverting Ha(z) gives the quantum corrected mirror map. We obtain for the first few 
orders in h and in Q expansion 

z{Q)W =z(Q)\ h s = 0, 

z(Q)W = -\Q 2 - 9Q 3 - 171Q 4 + . . . , (7.20) 

Let us now solve the difference equation via the WKB Ansatz given by (16. 3ft and 
(16.41) . The leading term of the solution reads 

f l-e x ±e~^e*(l-e*y-4 z' 
o {x) — log I - 

Higher order solutions S' n can be obtained easily, but are too length to be explicitly 
displayed here. Similar as in the previous example, expansion for z <C 1 effectively 
zooms into a C 3 patch of the geometry and the classical A-period can be read of from 
a contour integral around one of the singular points (c/., figure [2]), i.e., 

i(iff-i*M)-£f*. 
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Similarly, the higher order can be obtained via integrating S' n and we can indeed 
reproduce ( 17.19p . Let us turn to the B-period. The integration contour is indicated in 
blue in figure |2j We split the integration contour into two parts 

Taking the limit z <C 1 of the parametrization of the curve given in ( 17.18P restricts to 
a C 3 patch which includes /?i (c/., figure EJ). 

Explicit evaluation shows that for integrating along j3\ 

where A — > oo, 5 — > and we kept only finite terms. Let us calculate 



For example, for even n we obtain the first quantum corrections 

.2) 11 87 2 3349 o 176005 4 

U (2) = c z z 2 z 3 z 4 + ... , 

01 24 8 12 24 

U) 127 577 9 71081 „ 758401 , 

ni: ; = z z 2 z 3 z 4 + ... , 

01 5760 160 320 72 
with c some constant. In contrast to local P 1 x P 1 discussed in the previous section, 
we do not have anymore a Z 2 symmetry of the integration contour /3 and hence we 
can no longer expect that the integration along Pi and (3 2 yields up to signs the same 
result. Thus, we have to explicitly calculate the integrals along For that, we have 
to change the parametrization of the curve (I7.18P such that we end up in the limit 
z 1 in the patch which includes @2- It is not hard to see that the transformation 
x — > —x,p -4 — p combined with a change of complex structure achieves this. We 
obtain the curve 

- 1 + ze~ x + e~ p + e h/2 e x+p = , (7.21) 
which we quantize and solved via a WKB Ansatz as above. We obtain 

Hence, adding up the integrations along /3\ and (3 2 yield indeed the classical B-period. 

Integrating along /3 2 for some higher S' n yields 

, 2) 5 57 , 2509 _ 141355 . 

= c z - —z 2 z 3 z 4 + ... , 



12 4 6 12 

, 4 ) 67 317 2 45401 532285 , 

Yl = z z z + 

2880 80 160 36 
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Adding up IT^ = IL^J + IlL" and inserting the quantum corrected mirror map (I7.20j) , 

(2) 

we finally obtain (up to the constant part of U B ) 

(4) _ 29 207 55341 ^ 
VYb ~ 640 W 32 W 160 W + '- - • 

Let us compare this result to what we know. From the results of [I] one expects that 
the refined B-model 1-loop amplitude is given by 

J (/?) = - log (r) + — log (z) — — log (1 - 272) , 

with r := QdQz(Q)\f i o and one can easily infer that this indeed matches the corre- 
sponding refined vertex results. Using the usual holomorphic anomaly equations, it 
is straightforward to calculate in the B-model higher genus amplitudes J-( 9 >((3), which 
can again be matched with refined vertex results [BE] (this also has been confirmed in 
[6j). Comparing with ("17. 22[) shows that 

= Qd Q \im(3 2 ^ 2 \(3). 

8 Brane Wave-Functions and Liouville Amplitudes 

In section [4] we have seen that the brane partition functions of the /3-ensemble ma- 
trix models satisfy an operator equation which takes the form of the conformal Ward 
identity equation. When specializing to the Penner type logarithmic potential, it is 
well-known that the matrix model computes the Liouville conformal blocks with back- 
ground charge Q = b + b^ 1 , where 

b 2 = -P = * . 

In this case, the t\- and e2-brane are exactly given by the two types of degenerate 
vertex operators, denoted by V_j_ and V_b respectively, and the differential equation 

26 2 

( 14. lip becomes nothing but the usual BPZ equation [15] for conformal blocks with 
degenerate field insertions 

(b ±2 dl - £ —d Za - Y: j^^jZJx) = , 
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where A Q = a(Q — a) is the conformal dimension of the operator with momentum a. 
We have also argued that, provided we know the answer for the open string partition 
function, this knowledge can help us to deduce the closed string partition function, 
especially in the NS limit. We will see below how this idea is realized in this case. 

The key observation is that the ei-brane partition function ( 14. 2 p is always a polyno- 
mial of degree N in the brane location x. To be explicit, let us choose the logarithmic 
matrix model potential —f-W(x) = 2b Yla=o a <* log (a; — z a ) corresponding to the n- 



point function (V^ n _ 1 (cxD) Ila=o Va a ( z a)), where the vertex inserted at infinity carries 
momentum given by the conservation rule Yla=o a a = Q — Nb. Suppose we are inter- 
ested in the corresponding n-point conformal block computed by an appropriate choice 
of contour following [69J, we can probe the n-point function by inserting an ex-brane, 
or V_±(x) in the Liouville language, and arrive at the brane partition function 



Z 1 {x) = H(z a - z b ) 2a ^ ]J(x - z a ) aa/b 

a<b a 

x I d N z J] {zi - z 3 )- 2h2 Y[( Zl - z a ) 2a « b Y[( Zl - x) (8.1) 

l<i<j<N i,a i 

which is explicitly a polynomial in x. 

As alluded above, this ei-brane has the property that we can move it close to any 
other point z a of insertions without introducing any non-analyticity, and this gives us 
various relations among the n-point function with momenta shifted in units of — ^. In 
this way we can obtain essential information about the n-point function (the closed 
partition function) from the knowledge of the n + 1-point function (the brane partition 
function). 

In more details, consider the n-point function 

/N n-1 
d N z n & - z 3 y 2b2 n n( Zi - Za ) 2 ^ h , 

a<.o l<i<j<N i=l a=0 

it is easy to see that it satisfies the following relations to Z\{x) given in (18.11) : 

Z x (z a ) = lim Zi(x) = Z N (a a - ±) , (8.2) 

and 



Z dosed {a) = Ufa ~ Zb? aaab j d N z H (zt - Zj y 2b2 H 

a<b l<i<j<N i,a 
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Thus we find a set of difference equations for log Z c i ose ^ one for each non-degenerate 
insertion point, which depends on the open partition function Z\{%) but not on its 
normalization 

Z dosed (a a - ^)/Z closed (a a ) = Z^Za) / Zx(oo) . 

On the left hand side, we keep everything fixed in both the numerator and the de- 
nominator, apart from the momentum of the insertion point Z\ (x) approaches! From 
these, at least in principle, we can determine (up to normalization) Z c \ osed from the 
knowledge of Z\[x). We will now illustrate this procedure with the of 3 + 1-point func- 
tion example. In this case the potential W (x) has only one critical point and hence the 
corresponding spectral curve has genus zero. Recall that from the general discussion 
in section [6], the time dependence of the Schrodinger equation is trivial and one does 
not need to take the NS limit. 

In this case, the open partition function 

Zi(x) = x A -o-i/26- A -o- A i/26(i - x) a "i- 1 / 26 - Aq i" Ai / 26 ^ w (x) 

with 

r N 

*n(x) = d N z n & - z 3 y 2b2 Y[( Zl ) 2aob ( Zl - l) 2a ° b Y[( Zl - x) 

l<i<j<N i=l i 

satisfies the differential equation 

2 d 2 + A + 1 \ _ ^ao _ A Q1 + A ao + A Ql + A-i/26 - Aq 2 \ = ( 

dx 2 x x — ldx x 2 (x — l) 2 x(x — 1) / 



b 

with 



a + a 1 + a 2 -l/2b = Q - Nb . (8.3) 

This is solved by 

q? N ( x ) = cF(A,B;C; x) 

with 

A = -N , B = N-1 + \ (« + ai -\),C= |(a - ^) , 
and F(A, B;C;x) is the hypergeometric series 

(A) n (B) n x n 



F(A,B;C;x) = Y, 



(C) n n\ 



9 From the point of view of Liouville conformal blocks, Zi(oo) is related to 
(K* n _i-i/2&(°°) Ila=o K> a ( z o)) where the momentum of the insertion at infinity has been shifted by 
-1/26. 
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with (A) n = A(A + 1) . . . (A + n — 1). The normalization constant c remains to be 
fixed. Since A = —N is a negative integer in our case, the solution ^n(x) is a finite 
polynomial in x of degree N, as we argued earlier. From this, we can find the recurrence 
equations satisfied by the chiral half of the Liouville 3-point function 

r N 

Z N (a , 



( d N z n - z 3 r 2b2 fi(^) 2aob - 1) 20 

l<i<j<N i=l 



Using ( 18. 2D and colliding x with 0, 1 and oo and using analyticity of \&jv(x), it imme- 
diately follows that 

^(qp + l/2Mi) = ^tv(O) 
Z N (a ,a 1 ) ipN (oo) 

Z N (a ,a 1 + 1/26) _ ip N (l) 
Z N (a ,ai) iPn(oo) ' 

Using the special values 



F(A,B;C;x = 1) 
lim x~ N F(A = -N,B;C;x 

x—>oo 

the above recursive relation reads 



T(C)T(C - A - B) 
T(C — A)Y(C — B) 



(C) 



N 



Z N (a - ±, a ± ) j-j 1 ( 2ba - 1 + jb 



Z N (a , at) jj^ \2ba 2 - 1 + jb 2 

Z N (a , ai - ^) f-lba-y + 1 - jb 2 

Z N (a ,a l ) f} Q V 2ba 2 -l+jb 2 

where a 2 is given by «o,i as in ( 18. 3p . 

From this we conclude that the chiral three-point function must satisfy 

7 , ^ ft r(l-2to 1 - J 6') 

Indeed, after fixing the value Zn(0,0) as the initial condition, we then obtain the 
expression 

7{ , jt T(-(l+j)b 2 )T(l~2ba 1 -jb 2 ) 
N{ °' lJ Al T(-b 2 )T(2ba + jb 2 )T(2ba 2 + j& 2 ) ' 

which indeed gives the chiral half of the Liouville three-point function [TOJ [69] . 
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9 Conclusion 



In this work we have discussed several aspects of refined topological strings, emphasiz- 
ing the role of the branes. In particular, we have argued that brane partition functions 
are quantum mechanical wave-functions satisfying multi-time dependent Schrodinger- 
like equations, with times given by non-renormalizable moduli. The derivation was 
done using a matrix model realization of the refined topological string on specific ge- 
ometries. However, we expect this to hold in more generalities, though in general the 
precise identification of the times is not immediately transparent yet. 

In the NS limit, the time dependence drops out and we end up with a time- 
independent Schrddinger equation, making contact with the earlier results of [HI H2J 03] . 
The refined partition function in this limit can be recovered from the brane wave- 
function via considering the monodromies around the cycles of the local geometry, 
which we have also checked explicitly in several non-trivial examples. Using this fact, 
we have explained the observation of Nekrasov and Shatashvili in connecting integrable 
systems to gauge theory partition functions: The integrable system arises in the study 
of the world-sheet amplitudes and their target interpretation can be phrased in terms 
of gauge theory partition functions, thus explaining the NS result. 

There are various extensions of the present work one can consider. The simplest 
one, which should be straightforward, is to generalize the discussions in this paper from 
a single matrix model to Toda-like matrix models. This should in principle allow us to 
get arbitrary B-model geometries with arbitrary analytic H(x,p), where the powers of 
p, x are bounded. 

A more important extension of this work involves a deeper understanding of the 
brane partition function away from the NS limit. As already mentioned above, even 
away from this limit we expect an interesting wave-function which is now time- 
dependent. Uncovering the meaning of this wave-function from the perspective of 
integrable systems could be very interesting. Also, finding a way to compute via the 
time-dependent wave-function the complete refined topological string partition func- 
tion, as we have done for the NS limit via the time-independent wave-function, would 
be clearly important. We are currently elaborating on these ideas and plan to report 
on them elsewhere. 
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